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ABSTRACT 


When two or more bodies are in contact the load 
transfer takes place through the contact areas which 
are very small usually and the stresses in the contact 
region are very high. Plastic flow as well as fatigue 
failures may occur due to these high amount of stresses. 
Boundary Element Method is particularly suitable for 
contact problems and has been used to analyse such 
problems, whose boundary conditions are non-linear. Both 
static and rolling contact problems are analysed. Linear 
elastic bodies in two dimension are analysed with plane 
stress or plane strain assumption. Constant elements and 
linear elements are used for the numerical solution. Friction 
in the contact region is taken into consideration. Boundary 
conditions for static and rolling contact problems are 
discussed, A method to determine the contribution of micro- 
slip in the contact region towards the overall rolling 
friction coefficient is discussed. 
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CHAPTER 1 
INTRODUCTION 

Contact problems and the study of the load transfer in 
mechanical assemblages like rail and railwheel, cam mechanism, 
journal bearings, etc, are of great importance due to the 
severity of the stresses in the vicinity of contact region. 

If the stress exceeds the linear range there will be plastic 
flow. Moreover if cyclic stresses are applied then to find out 
the fatigue strength of the material, an analysis of the contact 
region becomes imperative. 

There are two types of contact problems: (i) Static 
contact problems and (ii) Rolling contact problems. Several 
analytical solutions for half planes are available. 

Hertz solutions [1 ] are for static conditions without friction, 
whereas solutions by Saiith and Liu [2] consider friction also 
in the contact region, Poritsky [3] considers friction and 
discusses the boundary conditions in the contact region for 
rolling contact, A review by Johnson [4] gives a comprehensive 
account of the nature of contact problems. 

Later many other solutions had been obtained using 
numerical methods like Finite Element Methods (FEM) [5,6,7 ] 
and Boundary Element Methods [8], In FEM most of the available 
solutions are for half-plane problems [5,7] , In all these 
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cases traction distributions along the contact region and 
transverse- to the contact region are found out by an iterative 
procedure in which starting traction along the contact region 
is chosen to be zero. In the subsequent steps it is increased 
slightly and the resulting traction distribution along the 
contact region is. supplied as the new distribution. Iterations 
are continued till tractions converge. Linear programming 
techniques were also used by treating the qontact conditions as 
constraints [ 6] , A solution of a contact problem is basically 
a solution of two or more Navier’ s equations with contact 
conditions. With appropriate boundary conditions both static 
and rolling contact types of problems can be solved. 

Boundary Element Method (BEM), a method which finds 
increasing applications in special problems like infinite 
domain, fracture mechanics, etc., can also be used for contact 
problems [8,9] , The method typically starts with weighted 
residual statement of the equilibrium equation. The weighting 
functions are Green’s solutions of the equilibriim equation. 
These weighting functions are of such nature that the domain 
terms in the we i^ ted residual statement disappears and only 
the boundary terms remain as an integral equation. Similar 
integral equations are obtained for each of the bodies involved 
in the contact. These equations can be solved using the contact 
conditions t© obtain the unknown values on the boundary. Using 
these boundary values internal stresses, if necessary, can be 
calculated. 
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BEM has the following distinct advantages over the 
domain methods such as FEM: 

(i) BEM treats only the boundaries which are of primary 
interest in the solution procedure, 

(ii) It is possible to couple normal and tangential tractions, 

(iii) The contact pressuresare obtained directly from the 
tractions which are primary unknown quantities and are 
determined with the same accuracy as the displacement 
unknowns, 

( iv) Sub-surface stresses are calculated, if they are needed, 
and they are also more accurate than displacement FEM 
formulation. 

In this work an attempt is made to determine the 
rolling friction while circular wheel rolls on a horizontal 
ground. The rolling friction is due to several mechanisms 
of energy dissipation such as micro-slip in the contact 
region, hysteresis losses, impact, etc, [ 10] , The present 
work attempts to determine the rolling friction due to 
micro-slip [ 1 1 , 1 2] in the contact region. Contact problems 
in 2-D with linear materials are analysed. Mainly cylinder 
on an elastic foundation is considered, with different 
combination of materials under different loading conditions. 

In Chapter 2, BEM formulation for isotropic, two 
dimensional problems has been described. Chapter 3 describes 
the BEM formulation for contact problems and presents the 
various sets of contact conditions. The solution procedure of 
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the contact problems is also discussed. Chapter 4 discusses 
the results while Chapter 5 gives conclusions. 
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CHAPTER 2 

BEM FORMULATION FOR ISOTROPIC BODIES 


The procedure of applying Boundary Element Method to 
linearly elastic problems is given in the following sections. 
Two dimensional isotropic bodies are considered unless other- 
wise stated. I^merical implementation and the method of stress 
evaluation are also presented. 


2,1 Fundamental Solutions ; 

Equilibrium equations for two dimensional isotropic 
bodies in the absence of body forces can be written as 

"jk.j = 0 (2.1) 


where 



the stress tensor. 


If 

conditions 


r is the boundary of the domain then the boundary 
, displacements u, and tractions p., are 


u j = u j on j = 1 , 2 


p. = p, on 3 = 1»2 

J J ^ 

(2.2) 


and r = Ej + T2 , 


Fundamental solution for Eq, 2,1 is the solution of the 
following equation: 


0 


(2,3) 



6 


where 6^ represents a unit load at point * i‘ in the 
direction 1 (Dirac delta function) and cr., is the stress 
tensor (Fig^ 2,1), So, the fundamental solutions for the 
two dimensional isotropic bodies (plane strain) are [13,14] 

^Ik “ ( 1- V) r ^Ik ^*1 ^'k^ ( 2.4) 


-’Ik = 4i( 


{-v)r ' ^an ^Ik 


-(1-20^) (r,j^ - r,j. n^) ] 


(2.5) 


where 


* 

Plk - 


'^Ik - 

G - 

V 

r - 



tractions in k directions due to a unit 

force in the 1 direction 

displacements in k direction die to a unit 

force in the 1 direction 

shear modulus 

Poisson’ s ratio 

distance between the point where unit load 
is applied (source point, i) and an 
arbitrary point in the field (sampling 
point, a) (Fig, 2,1) 
ccaaponent of r in k direction 
unit normal on the boundary (Fig^ 2,1) 
normal in k direction on the boundary 
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Physically these fundamental solutions are the 
resulting displacements and tractions in k direction due 
to a unit load applied at a point i in an infinite body 
in 1 direction [ 15 ] , 


2,2 Weighted Residual Statement : 

Considering the unit force is applied in one particular 
direction, the principle of virtual displacements for linear 
elastic problems can be written as [13 ] 


/ 

Q 


^3K3 



4 tPic - pp “k ^ 


( 2 . 6 ) 


» 

where u^ are the virtual displacements identically satisfying 

* 

the homogeneous boundary conditions Uj^ = 0 on r 

If we now interpret Uj^ as weighting functions which 
do not satisfy these conditions on then the expression 
can be rewritten as, 


/ 

Q 





( Pk-Pk) 


“k 


dr + 


/XUj^-Uk)Pk 


dr 

(2.7) 


where 





2,3 Boundary Element Equation ; 

Integrating Eqn, 2,7 by parts twice and substituting 


Eqn, 2.3, we get 





"l + f "k Pk Pk “k <Jr (2.8) 

where Uj represents the displacement at point i in the 1 
direction. 

If we consider unit forces are acting in the two 
directions (three directions for 3-dimensional problems) 
the above equation becomes 

"i + ^I,“k Plk '‘r=/PkPlk (2.9) 

This equation is the basic relation for developing the BEM 
models. In this equation the first term provides displacement 
of a point in the domain whereas all other terms relate to 
the points on the boundary, Eqn, 2,9 can be modified so that 
the first term refers to the displacement of a point on the 
boundary by assuming a small hemisphere around the point and 
letting the size of the sphere go to zero [3] , Applying 
the limiting procedure Eqn, 2,9 becomes 

(f + / Uj^ p*j, d r = / pj^ u*j^ dr (2.10) 

where the value of cf depends on the smoothness of the 
boundary, ^ contributes to Eqn, 2,10 only if the source 
point ' i’ and sampling point ’a*. Fig, 2,1, coincides, 

2,4 Numerical Implementation : 

The boundary is discretised into certain number of 
elements which are one- dimensional lines in the case of 2-D 
regions ( surfaces in the case of 3-D regions) . The elements 
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maybe of the following types (Fig* 2,3)t 

(a) constant elements 

(b) linear elements 

(c) quadratic and hi^er order elements* 

For constant elements the variables (u, p) will ranain 
constant within each element and they will be referred with 
respect to the node which is at the centre of the element. 

For linear elements variables vary linearly and so on* 

Eqn* 2,10 in discretised form is 



N 

S 

j=1 


pj pik “k 



IN ^ 

3=1 li^ ^ 


1 = 1,2 
k = 1,2 


where N is the number of elements, 

r. element ’j' on the discretised boundary (Fig*2.2), The 

%J 

summation in the second and third terms refer to summation 
over all the elements on the boundary, Eqn* 2,11 is a set 
of 2-squations per node, corresponding to each of the 
degrees of freedom. In Eqn* 2*11 each terra is a 2x2 matrix 
and 

First term = 

where [c^] can have values only if i =3 and the values vary 
according to the smoothness of the boundary, for example, 
on a smooth surface c^^ = = 0,5 and c ^2 “ ^21 




i) Constant 




.2.3 Different types of elements. 
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Second term (of Eqn. 2,11) for one particular 


f 

P11 

* 

Pi 2 




* 

P21 

* 

P 22 

< 

f 

^2 




Third term of Eqn, 2,11 for one particular j 



u 


u 


11 

21 


u 


u 


12 

22 


P2 


d r 


For constant elements U|^ and pj^ ( k = 1,2) will remain 
constant in an element. In case of linear elements Uj^ 
and pjj, will be varying within the element. They can be 
expressed using the nodal values and the shape functions as 


Pk 


0 Uj^ 

T t 

^ Pk 



(2.12) 


where. 


I 



nodal displacements of the element 



nodal tractions of the element 

(0^ and 02 shape 

functions for linear 
element) 


01 02 0 0 
0 0 0 | 02 
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"til 

Referring to Fig, 2.1, and q. are the nodes of j-1 

element and q. and q. . are the nodes of the elaaent. It 

J J + 1 

can be observed from Eqn, 2,11 and Eq, 2,12 that the variables 
connected with q. will figure in the integration of both 

+K 

j -1 and j elements. Similarly all the nodes will 
contribute to the adjacent elonents. So, for each i node 
when the integration is carried over the boundary a set of 
two equations with 2N variables under each integral term are 
obtained. For each node Eq, 2,11 gives a 2x2 matrix associated 
with the variables of node under each integral term. In 
matrix form Eq, 2,11 for source point i is 
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PI “ 


are displacements and tractions of node i. 


It can be noted that c term has been added with the 
value given by the integral term in LHS when i = j (Eqn, 2,11), 
The evaluation of c^ will be explained later. The submatrices 


associated w 


ith Jl!l I and J^ll - both are the variables of 

1^2/1 tP 2 jl 


node 1, are [h] = [9] = 

"21 22 


^11 ^12 
^21 ^22 


Similarly all the nodes have such sub-matrices (Eqn, 2,13), Wien 
Eqn, 2,11 is evaluated for all the boundary nodes in turn, 
matrix equation shown below is obtained: “ 


[H] lu } = [G] Cp} 


(2.14) 


where, 

[ H] — a matrix of size 2Nx2N with N-subma trices [ h] 
of size 2x2 

[ G ] — a matrix of size 2Nx2N with N-submatrices [ g ] 

of size 2x2 

{u land {pi- displacement and traction vectors as given in 
Eqn. 2,13, 



Inspection of fundamental solutions, Eqn, 2,4 and 
Eqn, 2,5, shows the singularity arises in Eq, 2,11 when r 
becomes zero, i.e,, v\iien source point ’i’ itself is also the 
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sampling point. These singularities (i = 3 ) contribute to 
the diagonal submatrices of size 2x2, Eqn, 2,13, to both [ H] 
and [G] in Eqn, 2,14, When i =3 the first term of Eqn, 2,11 also 
contributes to the diagonal submatrices of [H] . Being a 
logarithmic singularity right hand side integral (.when i = 3 *) 
will be analytically evaluated in the principal value sense 
without much difficulty. In case of diagonal submatrices of [H] 
no analytical method is attempted, but they are evaluated from 
rigid body motion principles. By giving unit displacements in 
all directions Eqn, 2,14 becomes 

[ H ] C U } = 0 (2,15) 


where CU} is a vector defining unit rigid displacement in all 
the directions. So, the diagonal submatrices of [H] will be. 


[hy] 


m 


S 

k=1,N 

lq4m 


t"i3^k 


(2,16) 


i= 1,2 

3 = 1,2 

where 

N - the number of nodes 

m - the node which gives rise to the singularity. 


All other non-diagonal elements of both the matrices are 
numerically evaluated using ’Gaussian Quadrature’ , Hence 
all the coefficients of [H] and [G] will be known. 
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It may be noted that in Eqn, 2,11 out of the 4 variables 
uj^ and pj^ ( k = 1,2) at each node two will be known boundary 
conditions. If there are N nodes after discretization there 
will be 2N unknovjns associated with them. Bringing all the 
known values to right hand side and the unknown values to the 
left side Eqn, 2,14 becomes 

[A ] C Xl = C B } (2,17) 

v\here 

[ A ] - matrix of size 2Nx2N 

{X } - unknown quantities of Uj^& Pj^ in the form of 

vector of size 2N 

C B } - column vector of size 2N, calculated from 

the known boundary values. 

Solving Eqn, 2,17 with any standard method the unkno\A?ns 
of the Uj^ and Pj^' s on boundary are obtained and thus all the 
boundary values are known, 

2,6 Determination of Stresses in the Domain ; 

After solving the system given by Eqn, 2,17 we will 
have all the boundary values, i,e,, displacements and tractions. 
With the help of Eqn, 2,8 we can get displacements at any point 
in the domain. Using these inter ral points, we can find out 
strains and hence the stresses. 
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Alternate way to determine stresses directly is to 
differentiate Eqn, 2,8 at the internal point. Substituting 
the fundamental solutions in the resulting expression we can 
get an expression which will directly give us the stresses as 


[13,14] 





P„ d r 



U, 


d r 


where 



(2,18) 


D 


kij 


“ 4in:( 1-.y)r ^ 1 - 20 ^) { 


- 6 


ij 


+ 2 r,i r,, r.i^ } 


®kij 


^ -5 r,^ + r,,) 


4x(1J-)r^ -- an - ■ - ' ij *'k -“ik *>j — Jk 

-4 r, ^ r, ^ r, j^] + 2v( n^ r, ^ k ” j i k 


+ (1--2y)(2 njj ^*i ^j^ik '^i ^jk^ 

- ( 1-4v) nj^ 6 } 


The above equations can be used to evaluate the stresses at 
points within the domain only. If the points are on the 
boundary a modified approach is adopted as explained in the 
Appendix A, 
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CHAPTER 3 

FORMULATION FOR CONTA:;! PROBLEMS 

In this chapter contact problem is defined in the first 
section. After that contac t conditions related to static 
contact, with friction and wiihout friction, are discussed. 
Contact conditions for rolling contact are presented. Solution 
procedures for both types of contact problems are given. Two 
dimensional problems with linear materials are considered, 

3,1 Definition of the Contact Problem; 

Consider two linearly elastic bodies A and B bounded 

by cylindrical surfaces. The traces of the boundary surfaces 

A B 

in the x-y plane are the curves, r and r respectively. 

The geometry and all variable quantities are independent of 
z-coordinate. The problem is thus considered as a plane problem 
in the x-y system Fig, 3,1, The tractions acting on the 
boundaries are denoted by t^^ and the displacements by Uj^, On 
the boundary r the two bodies are in contact with each other. 
The parts of the boundaries corresponding to the contact 

A 

boundaries r are in the undefomed state denoted by r ^ 
c c 

D 

and r“. ■’’’^2 natural co-ordinate system with 

along the normal and TI 2 along the tangent of the point on 


the boundary 
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3,2 Contact Conditions ; 

Broadly there are two types of contact problems, viz, 
static and rolling contact. Contact conditions for both 
types are given below. 


3,2,1 Contact Conditions for Static Contact : 

The contact conditions are introduced in natural 
co-ordinates and later all the equations will be converted to 
global co-ordinates. 

Displacements u-j and ti 2 belong to the undeformed 
surfaces in natural co-ordinates (Fig, 3,2) whereas v-j and 
are displacements of the same surfaces after contact has been 
established, 1 and 2 represent the directions of natural 
co-ordinates. Relations between these two displacements are 
given below [s] ; 




- a on 

- (1-a) u° on r® 

A 

on 

on 


(3.1) 


where a - a function of the location of contact boundary 

u°- the distance in direction 1 between the original 
positions of the points that are now in contact. 
Superscripts A and B indicate the respective bodies. 
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Static contact conditions can be given for three 
situations viz,, 

(i) No-friction 

(ii) Infinite friction 

(iii) Finite friction 


(i) Contact conditions for no-friction case can be stated as 
follows on 


V 


1 


+ V 


e - 


B 

1 

B 

1 

A 

2 

B 

2 


t^ ^ 0 , 


= 0 

= 0 
= 0 
= G 
t^ < 0 


(3.2) 


( ii) For infinite friction on r 


v^ + 


V| = 0 

+ Vrt = 0 

< - t? = 0 

tg - ll = 0 


(3.3) 


(iii) For finite friction 


Tea 


V., + 


B 


.B 


h” - - 

S - 

4 “'^^2 

4 - 4 = 


0 

0 

0 

0 


if I "^2 ^ ^ ^ 


(3,4) 
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Tcs = *2 = ± n t]' k = A.B 

.A .B 

t2 - 0 

if 1 12 I ^ p. It^ I 

where 

- contact zones 
“ stick zone (No-slip) 
r^Q - slip zone 


Sign of p, coefficient of friction, is decided in such a 
manner that 

sign (t 2 ) sign (v^ + v®) 


The idea is that p should have that sign which resulted in 
energy dissipation from the earlier configuration. 


Using liie relationship in the form of Eqn, 3,1, 


A B 

vi + 


i,e. 


Uh 


A ^ . B 
B 


= u^ + u^ - 


= G 


+ u 


1 


= u. 


Similarly, 



B 

+ V2 



If we assume the contact length or area is very small 
and if it is almost parallel to the x-axis, then Eqn, 3, 2-3, 4 
can be rewritten in the global co-ordinates in the follovdng 


manner: 
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(i) No-friction; 

0 
0 
0 

(ii) Infinite friction: 




(iii) Finite friction; 

A 


ca 


GS 


u. 


“y " 

uA 




B 

tt = 
X 


u 

1 

0 

0 

0 


if I t^l < P Ityl 


% = 


X 

t^ + t® 

X X 


y 

= 0 


if itjji flirty I 


(3,5) 


(3.6) 


( 3.7) 


Subscripts x and y denote the directions x and y of global 
co-ordinates. 
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In this work contact conditions given by Eqn,3,5 - 3,7 
are used. 

If the contact surface is sufficiently large then all 
the conditions given in natural co-ordinates are converted to 
global ones using transformation matrices [8] , On the other 
hand, the ehtiire problem can be solved in natural co-ordinates 
itself [ 8 ] , 


3,2,2 Contact Conditions for Rolling Contact ; 

Unlike in static contact problems, no friction and 
infinite friction cases are not considered. Only finite 
friction case is dealt with as it is more practical than 
the other two conditions. 


Analytical approaches [3,12] based on half-plane solutions 
suggest that the leading edge of the contact zone should have 
slip region. However experimental studies [16] indicate 
that there will be indeed slip zone on both the sides with 
no-slip zone in between. So, the contact conditions for rolling 
contact are as given below: 
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where is the term that takes care of the rolling contact 
situation. This ’relative displacement', 6 , is found out 
by an iterative procedure, wAiich will be explained in detail 
later in Sec. 3,4, 

3.3 Solution Method of Contact Problems by BEM 

Consider two linear elastic bodies A and B in contact 
(Fig, 3,1). The boundary element equation for these two 
bodies is 

“i + ri = / «ik Pk " (3-9) 

m = A,B 

In ordinary stress analysis problems half of the 
boundary values will be known and the renainings are found out 
by solving Eqn, 3,9 as explained in the previous chapter. In 
case of contact problems, neither displacements nor tractions 
will be known in the contact region and so the number of 
equations available will be less than the number of unknowns. 
Sc, we have to introduce certain equations which govern the 
contact phenomena and such equations are the contact conditions 
discussed in the previous section. In this manner number of 
unknowns and number of equations available can be made equal 
and the probloa is made solvable. 

Initially integration of Eqn, 3.9 is carried over the 
boundaries of A and B independently and the resultant matrix 



equation is (refer Eqn, 2.14) 
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bodies A and B respectively. 


Variables of contact conditions are rearranged such 
that in each equation variables of the upper body (A) are 
expressed in terms of that of lower body (B). According to 
these relations, in Eqn. 3.10 columns corresponding to the 
degree of freedom of the nodes of the lower body are added to 
the corresponding columns of the upper body. By appropriate 
matrix operations all the unknovsn variables are brought from 
the right hand side to the left hand side matrix. After this 
column transfers for the applied boundary conditions are 
carried over as explained in Sec, 2.5, The right hand side 
column {t } has only known values and therefore Eqn. 3.10 
acquires the form [A] CXI = {B} which can be solved by any 

standard method available. 


3 ,4 Iteration Procedure to Determine Contact Zones for 
Pinlie P'riciion ; 

The contact conditions are of non-linear as the contact 
length will not grow linearly with the applied load so, as 
any other non-linear problem, to determine contact length 
iterations have to be done. Iteration procedure for static 
and rolling contact problens are separately given below. 
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3* 4, 1 Static Contact Problems : 

To start with certain portion of the boundary is 
considered to be in contact and all the necessary contact 
conditions are implemented. The resulting system of matrices 
is solved to get the displacements and tractions throughout 
the boundaries of both the bodies. Since in the contact 
region there could only be the compressive stress in the 
direction transverse to the boundaries of contact, this is 
taken as a reference to decide the contact region, Fig, 3, 3, 

In all the assumed contact nodes, there should be compressive 
stress, otherwise those nodes which have tensile stress will 
be rauoved from the contact region. If all the nodes have 
compressive stress, few more nodes are added to the contact 
region and the corresponding contact conditions are introduced 
and the resulting system is solved. This process is repeated 
till the extreme nodes just begin to have tensile stress. 

Then with slight rearrangement of the extrane nodes, the points 
where there are neither compressive nor tensile stress are 
located. This contact region is the region of contact for 
the given load and materials. Refer Fig, 3,3, 

In determining the ’zero pressure' points, different 
criteria are followed for constant and linear elements. 

In case of constant element, the extreme pair of 
contact elements of the contact region on both the bodies 
are reduced in the length as soon as extreme elements become 
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Too large contact area 
Positive normal stress 



Correct contact solution 


4 



Too small contact area 
Geometrical incompatibility 


Fig. 3.3 Three possibilities of solution as results 
of different choice of contact area. 



Fig. 3.4 Relative position of the nodes in the 

contact region in rolling contact problems. 

(Gap magnified for clarity) 
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tensile, A better approximation to the contact length is 
obtained by interpolating using^the tractions in y-direction 
of the extreme pair of elements. With this rearrangement 
entire solution procedure is carried out again. This procedure 
is repeated with the updated contact length till a satisfactory 
value for contact length is obtained. 

In case of linear elements, when the extreme nodes 
acquire tensile stress, interpolation is carried over between 
extreme pairs of nodes and on both the bodies the lengths of 
extreme elements alone (unlike in the previous case #iere a 
pair of extreme elements together were adjusted) are adjusted, 
i,e,, two elements on the upper body and two elements in the 
lower body and hence 4 elements are adjusted as against 8 
elements in the case of constant elements. 

If finite friction is to be considered in the contact 
region, to start with it is assumed that there is no slip in 
the contact region at all. After establishing the contact 
region in the above manner, nodes in which the condition 
lt„l xplt I is not satisfied are located. For all these 

X y 

nodes contact conditions corresponding to slip- region 
( Eqn, 3,7) are introduced. With these new set of contact 
conditions matrix system given by Eqn, 3,10 is solved again. 

Once more all the contact nodes are checked for the violation 
of the condition I tj^I < l and if there are such nodes, 
they are joined with the slipping nodes category and the problem 
is solved again. This is continued till all the nodes obey 
the condition * < V- I U 
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After that, the nodes immediately outside the 
contact region are examined for interference, i.e, , one 
body may be piercing the other one. In such case even 
those nodes are added to the contact region and the problem 
is solved. On the other hand if there is no interference, 
the contact length already found out is taken as the correct 
one and the corresponding traction distributions and displace- 
ment fields are used to find out the stress fields, 

3,4,2 Rolling Contact Problems ; 

All the techniques and conditions used for establishing 
the contact length in case of static contact problems are 
applicable for rolling contact problems also. However, there 
is an additional iterative loop vihich is explained in the 
following. 

Initially the probl^ is solved as a static contact 
one. While developing the slip zone in the contact region 
a relative displacement corresponding the nearest slip node 
(vAiich is in the leading edge of the contact length) to 
the no-slip region is introduced in the no-slip region also. 

In Fig, 3,4 thickened portions are slip regions while the 
region in between th®n is the no-slip one. For rolling 
contact problems the relative displacement corresponding to 
the node 4 is taken as 6 ^, Eqn. 3,8, if 1-2-3-4 is the leading 
edge. Using Eqn, 3,8 with this 5 ^ contact problem is 
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solved. This will give a new relative displacement for 
node 4, New 5 ^^ is taken to be equal to the nevj relative 
displacement for node 4 and the problem is solved again. 

This iteration procedure continues till the traction distri- 
butions in the contact region converge. If there is scope for 
adding more elements into slip region, iteration continues with 
new set of slip nodes and again the process explained above 
is repeated. So, determining § is an iterative loop inside 

A 

a general iterative loop to determine the slip region of the 
contact length^ 
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CHAPTER 4 

RESULTS AND DISCUSSION 


In this chapter results for thick cylinder v;ith internal 
pressure, static contact and rolling contact are given along 
with their discussion. All the materials are assumed to be 
isotropic and linear elastic. 


4,1 Thick Cylinder with Internal Pressure ; 

This is an axisymmetric, two dimensional, plane strain 
problem. Using symmetry, only one quarter of the thick 
cylinder is analysed. The boundary is discretised into 52 
elements and the elements are of constant type. Analytical 
solution for this problem is available for the stresses and 
the displacements in the body. Stresses obtained by BEM and 
the values of stresses obtained from analytical solutions are 
given in Table 4,1, The details of the problem are given below: 


Inner radius of the cylinder 
Outer radius of the cylinder 
Internal pressure 
Shear modulus 
Poisson ratio 


= 3,3 units 

=r 7,5 units 
=r 100 units 
= 94500 units 

= 0.1 


As can be seen results are in good agreement thus indicating 
the correctness of the code developed. The discretization 
is very fine and in fact a model with 26 elements also gave 
good results. 



Table 4,1 ; Results for Thick Cylinder 
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4,2 Static Contact Problems ; 

The BEM code is applied to solve static contact 
problems, A few problems solved by Ma [9] are solved for 
comparison purposes. They are semi-cylinder over elastic 
foundation and semi-cylinder over semi-cylinder. Loading 
conditions are as shown in Fig, 4,1, Constant elements are 
used. No friction case was solved in [ 9 ] , However, the 
results of this work are for finite friction. For these 
problems Hertz solutions are available for no-friction state 
in the contact region. The distribution of contact pressure 
is shown in Table 4,2, Agreement is satisfactory between 
the present results and that are given in [ 9 ] • If friction 
in the contact region is considered, then there will be 
transverse traction in the region. This region can be sub- 
divided into two regions viz, (i) no-slip zone and (ii) slip 
zone. In no-slip zone pt^ will be less than p t^ and pt^ 
will be equal to pt^ in slip zone, where t^ and ty are tractions 
in the contact region in x and y directions and p is the sliding 
coefficient of friction. Though there was not much change in 
the contact length with respect to p , lengths of slip and no- 
slip regions depend on p. Lower the coefficient of friction 
larger the slip region for a particular load and vice versa, 

A limiting value of p can always be found out such that the 
entire contact region may be having slip region alone, thus 
resulting in macro-motion of the body itself. 



Py = -40 N/mm^ 

i_LLi 

R =40nnnn 
))=0.3 

G = 83000 N/ mm^ 


-►X 


7 }i)j Di/nn )iniii 


Py= -10N/mm^ 



Fig. 4.1 Contact between semi -cylinder and 
elastic foundation and semi -cylinder 
and semi-cylinder. 



2) Inclined load passing through the origin 

Fig. 4.2 Cylinder on an elastic foundation. 
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Table 4 2(a) : Results of Semi- cylinder over Elastic 

Foundation (Ref, Fig, 4,1) 


2 

Contaot Pressure ( N/inm ) 


x(mm) — — r 

By Hertz 02] 


1,209 

0,0 

1,073 

700,0 

0,875 

1156,5 

0,625 

1470,1 

0,375 

1646,4 

0,125 

1727.7 


By Ma [12] 

By BEM 

0.64 

-0,3114 

705,6 

673,56 

1140,4 

11 34,3 

1484,2 

1460,4 

1657.8 

1631.7 

1739,3 

1713,6 
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Table 4.2(b) : Results of Semi-cylinder over Semi-cylinder 

(Ref, Fig, 4,1) 


2 

Contact Pressure (N/ram ) 


X 

(mm) 

By Hertz [[1 2] 

ByMa[12] 

By BEM 

1.192 

0.0 

0.9 

6,0 

1.621 

127.0 

128.9 

137.56 

1.283 

198.5 

200.3 

195,47 

0,977 

236.1 

243.6 

239,41 

0,672 

261,3 

259,1 

259,96 

0,367 

274.1 

279.1 

277,21 

0,122 

279.3 

281.8 

282,34 
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teMlts with Vertic al Load on the Cylinder on Foundation ; 

A long cylinder over an elastic foundation problem with 
vertical load (Fig, 4,2) was analysed. The width and height of 
the foundation are around 4 times and 3 times the radius of 
the cylinder respectively. Hertz solution is available for 
no-friction case and Smith and Liu's [2] solution is available 
for friction cases. 

Numerical results by applying BEM technique with constant 
elements is available in [8] , Fig, 4,3 gives the solutions 
from the present work and Ref, [8] , It may be noted that t^ 
distribution of both works are in agreoDoent whereas there is 
slight discrepancy in the case of t^^ distribution. In Ref, [ 8 ] 
loading is done in several increments whereas in this work 
total load is applied in one step because of the high amount 
of computer time involved. Nevertheless the traction distri- 
butions thus obtained are sufficiently accurate to give a 
satisfactory stresses as will be demonstrated in the later part 
of this section. 

Distributions of contact pressure ( ty) and t^^ in the 
contact region for a cylinder of radius 72 mm is shown in 
Fig, 4,4, Linear elements are used for modelling. This 
class of problems does not show much difference between friction 
and no-friction cases, regarding ty distribution whereas t^ 
distribution changes with coefficient of friction, because of 
the change in slip and no-slip zones. 



Vertical load =-590N/mm 
E s 4000 N/mm^ 
D = 0.35 
R = 50 mm 


p = 0.05 


□ - ty 

O- tx/p 
X - ty , 
A- tx/p 


Andersson [8 ] 


Present Work 


-5.0 -2.5 0 


-100 


5.0 . 


Contact Length, mm 


Fig. 4.3 Cylinder on elastic foundation, ty and 
tx for cylinder in the contact zone. 



- Contact pressure 

- Horizontal traction 
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For a vertical load on the cylinder. Fig, 4,2 , 
variation of stresses along the vertical diameter is shown 
in Fig. 4.5. Along with this plot solution by Smith and 

[^] also shown. While the solutions are matching for 
o'y and cr^y around the contact region, marked variation is 
seen in the case of cr^ near the contact region itself. For 
some points [2] gives as compressive stress whereas BEM 
solutions for these points give tensile stress values. This 
is due to the fact that Smith and Liu [2] solution was obtained 
from a basic half-plane solution. The BEM solutions can be 
checked in the following manner. Away from the contact region 
this solution will be similar to a cylinder loaded at dia- 
metrically opposite points [1 ] , Assuming each of the loads 
produces a simple radial stress distrioution, it can be found 
out what forces should be applied at the circumference of 
the disk in order to obtain such a stress distribution. Stresses 
on the boundary should be nullified by applying equal and 
opposite stresses, if the boundary is traction free. So, away 
from the point loads this problem and the contact problem are 
one and the same according to St, Venant’s principle. As 
can be seen from Fig. 4,5 away from the contact region BEM 
solutions and Hertz solution in[ 1 ] are matching well, 

A similar superposition procedure was carried out with 
elliptically varying load being applied instead of the point 
loads on the cylinder. Half-plane solutions for this kind of 
load is available in [2] . However, unlike in the point load 



of the solutions from different sources for a cylinder 
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case, here the stress field to be applied on the cylinder 

to make it traction free on the boundary is not uniform. 

So, solution for this stress field is obtained from BBS itself 

and used for superposition. The final solution of a and a 

y 

were in agreement with Smith and Liu near the contact region 
and with Hertz [1 ] away from the contact region. The final 
solution of O’ ^ was in agreement with BEM solutions at all 
the points. However, the results are not presented as this 
problem is almost similar to the point load case, 

4,2,2 Results of Inclined Load on the Cylinder ; 

When compared with the results of vertical load, there 
is marked change in the t distribution ( Fig, 4,6), This is because 
the inclined load that introduces a horizon’tal force and contact 
region is the only area where the reaction for this horizontal 
force will exist. In this case t^ distribution is also changed. 
It is skewed in such a manner that the resultant vertical load 
can counteract the moment introduced by the horizon’tal load. 

Stress distribution curves for the slightly inclined 
load (Fig, 4,2) are shown in Fig, 4.7. It may be noted that 
this load does not introduce any moment. Stress distributions 
due to vertical load 150,6 N/mm are also shown in the same 
figure for the sake of comparison. 




I 


Fig. 4.6 tx and ty distributions on the contact region of a cylinder 
loaded with an inclined load. 


Vertical load =-150,61 N/mm 
Inclined load = 150.61 N/mm 
Radius (R) =72 mm 

E = 4000 MPa 
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contact for vertical and inclined loads. 
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Contact or Rolling Contact Problpma ? 

A static contact problem itself is a highly nonlinear 
one and many such problems have to be solved with short time 
intervals to simulate rolling contact. So, for dynamic 
problems calculations involved are very complex. Computational 
work is reduced very much if only slow rolling is considered. 

It is assumed that contact conditions, region etc, are not 
affected by the value of rolling speed. Contact conditions as 
discussed in Section 3,2^2 are used. 

In order to study the effect of free rolling of a cylinder 

on an elastic foundation, loads are applied as shown in Fig, 4,8, 

For each normal load, various values of horizontal loads are 

applied. Horizontal loads are of small magnitudes with 

approximately 1/5GOO to 1/100 of the vertical load P and P 

X y 

correspond to the applied vertical and horizontal tractions 
whereas t^^ and ty correspond to the traction distributions at 
the contact region, 

4,3,1 Results for the Case of Same Material for Cylinder and 

Founclaiion; ^ ' 

When same material is used for both cylinder and foun- 
dation the distribution of ty and t^^ are as shown in Fig, 4,9 
for one particular applied Py and P^^, It can be noted that ty 
distribution in the contact region is skewed to the left. This 
is expected because of the moment introduced by the horizontal 



Z ‘151 
o 


V-tx/u 


Applied load =-32.8AN/mm 
Py =-6,0N/mm2 
Px =-0.005 N/ mm ^ 
E = AOOON/mm^ 
R = 31.4 mm 

M = 0.01 

V X) = 0.35 

V X-axis 1:0.1256 
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Contact Length, mm 
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load. The skewed distribution of t^ introduces a counter 
moment which can be calculated by using the t^ distribution 
itself. 

Stresses for points on the vertical axis, Fig, 4.10, 
are given for two horizontal loads and for one vertical load. 

Due to the asymmetry in the stress field wtfien compared 

with the solutions for vertical load alone, has lower values 
for all the corresponding points. As increases the point 
where maximum oy occurs also keeps moving to the left. For 
the same reason which was having zero values in the 

case of vertical load only along the line of symmetry, increases 
for the corresponding points, 

cr increases with increase in the applied horizontal 
load because of the nature of cr^ stress field (Fig, 4,11), 

On any given horizontal plane the trend of distribution 
is as shown in Fig, 4,11, As shown in the graph, this distri- 
bution when skews to the left value of at the centre-line 
increases with applied horizontal load, (Though the graphs in 
Fig, 4,11 are for Teflon cylinder and steel foundation, same 
trend is noted in the present case also), 

4,3,2 Results for the Case of Different Materials for Cylinder 

and Foundation ; 

The above study was also conducted for two other combina- 
tions viz, (i) PMMA (plexiglass) cylinder on steel foundation 
and (ii) Teflon cylinder on steel foundation. 






Distance in a Horizontal Plane, mm 


Fig. 4.11 Ox distribution on a plane y=2mm of a 

teflon cylinder rolling on steel foundation. 
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It is noteworthy of attention to see the difference in 
the pattern of t^^ distribution. Fig, 4,12, between this case 
and the former one. Fig, 4,9, where both the bodies are made of 
same material, A parametric study was carried out to see wJien 
this distribution changes its pattern from the one shown in 
Fig, 4,9 to the present one. The modulus of elasticity for the 
cylinder was kept constant and that of foundation was increased 
gradually. At a ratio of around 1,1 almost all the t^^ in the 
contact region vanishes and afterwards the present trend sets in. 
Actually the source of t^ is the interaction between displacement 
fields of the cylinder and the foundation. As the elasticity 
modulus of the foundation is raised there is more deformation 
in the cylinder and hence the change in traction distribution 
pattern. However, t^ distribution ( Fig, 4,13) retains the same 
trend as that for cylinder and foundation made of same material 
because it is mostly governed by the vertical load. 

Stresses are calculated along the vertical axis of the 
cylinder for both the combinations for a few vertical loads 
with different horizontal loads, a y. and stresses are 
plotted in Fig, 4,14 with vertical load as 32,84 B|/mm and 
horizontal loads as 0,0547 N/mm and 0,02737 N/nia for Teflon 
and steel combination. Results of a few other cases are given 
in Table to 4,3 to 4,7, 

In order to determine the contribution of micro-slip 
towards the overall rolling friction coefficient, we have to 
find the horizontal load necessary for a steady rolling. The 
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Fig. 4.12 Traction along x direction in the contact 
region for various horizontal loads acting 
on a PMMA cylinder rolling on a steel base. 
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Py= -10 N/mm^ 



Fig. 4.13 Traction along y -direction for various hori- 
zontal loads acting on a PMMA cylinder 
rolling on a steel foundation. 


Distance Along Vertical Axis, mm 


Applied Py=-6.0N/mm^ 

(Load=-32.84 N/mm) 


Compressive Stress, MPa 


Fig. 4.14 Comparison of stress fields for two 

different horizontal loads in a Teflon 
cylinder on steel foundation. 
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resultant moment about the centre of the cylinder due to 

r- 

each loading condition is calculated and for a particular, 
^-Vertical load when value of horizontal load is small this 
resultant moment will act such that a free rolling cylinder 
will decelerate , As the horizontal load is increased 
resultant moment gradually changes its sign so that the body 
starts accelerating [16], When the resultant mcraent is zero, 
that situation corresponds to steady rolling. For one parti- 
cular case of PMMA and steel combination, resultant moment is 
plotted against the applied horizontal load in Fig.4,15, 
Resultant moment is calculated by finding the horizontal and 
vertical reactions in the contact region and taking moments 
about the centre of the cylinder. 

The position of the resultant vertical force in the 
contact region shows very smooth behaviour, Fig, 4,15, with 
applied horizontal loads for all the cases, whereas the 
resultant moment is not having a uniform trend. These trends 
can be explained as follows. From Fig, 4,12 it can be seen 
that horizontal reaction will be the algebraic addition of 
positive and negative areas of a t^^ distribution curve. 

While a slight different distribution may also be accepted 
as a correct one, horizontal force in the contact region may 
vary considerably. In the process error gets accumulated 
in horizontal reaction. So, being a small value, it is very 
likely that the value of resultant moment will have large 
error and this is suspected to be the reason for an unstable 
behaviour of resultant mcmaent^ 
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vertical force in the contact region against applied horizonta 
load for PMMA cylinder on steel foundation. 
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From Fig, 4,15 rolling friction coefficient, p g, 
can be calculated if the horizontal load for zero resultant 
moment is divided by the applied vertical load and for PMMA 
cylinder on steel foundation \i^ = 0,0033 from the BEM results. 
This value is much higher when the usual value is around 
0,001 , Improvement in p can be expected if the numerical 
instability in the calculation of horizontal reaction force at 
the contact region can be avoided. 




Table 4,3 : PMMA Cylinder on Steel Foundation (Plane Stress) s Vertical lead = 75,3 N/ram. Traction 

Distributions on the Cylinder 

P„ s -6 Radius = 72 mm, 0,4, y= 0,35, = 4,57 GPa 



60 


Table 4,4 : PMMA Cylinder on Steel Foundation (Plane stress)* 

Vertical Load = 150*61 N/rnm, Traction Distributions 
on the Foundation 

Py = -12 Radius = 72 rm, \i= 0,4, v= 0,35, %>MMA = H 


x-co-ordinates 

mm 


= -0.01 N/mm^ 

t 

X 

t 

... . - 

1.692 

-0.73452 

-1.8363 

1 ,440 

-8*5447 

-21.362 

1,152 

-14.757 

-36*892 

0*864 

-18*414 

-46.036 

0.576 

-17.986 

-53*941 

0*288 

- 5,2529 

-58.538 

-0,288 

3.6150 

-59.217 

-0.576 

11,819 

-56,535 

-0,864 

20.155 

-50*382 

-1.152 

16,517 

-41 ,293 

-1.44 

15.245 

-33.111 

-1.8 

-0.031327 

- 0,078318 


Continued 
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Table 4,4 (Continued^ ; 


■coordinates 

mm 

Px = 

-0,025 N/mm^ 


"v 

1.566 

-0,52807 

-1.3202 

1.2528 

-12,519 

-31,297 

0,9396 

-16,505 

-41.263 

0.6264 

-20,375 

-50,937 

0.3132 

-8,8185 

-57.582 

-0,3132 

1 .9841 

-59,647 

-0.6264 

10,256 

-57,095 

-0.9396 

20.296 

-50.741 

-1 ,2528 

16,342 

-40,854 

-1,566 

12.34 

-30,850 

-1.8792 

-0.05596 

0,1399 


• • 9 • 9 • # 


Continued 


Table 4,4 ( Continued'^ i 


x-coordinates 

mm 


= -0.05 N/rom^ 


^x 

t 

y 


1 .4472 

0.61635 

1 .5409 


1.152 

-11.767 

-29.417 


0.8928 

-15.248 

-38.12 


0.576 

-19.264 

-48.159 


0.288 

-14.099 

-55,724 


-0.288 

- 0.59893 

-59,974 


-0.576 

4.9905 

-58.833 


-0.864 

13.647 

-55,7 


-1.152 

19.553 

-48.882 


-1.44 

15.817 

-39.543 


-1.728 

12.052 

-30.131 


-2.0304 

- 0.30543 

+ 0,76358 



X- coordinates 

= 

-0,1 N/EBn^ 





1,152 

- 0,37684 

-0.94211 


0.576 

-17.911 

-44.777 


-0,576 

- 4,9532 

-60,322 


-1.152 

1 3.727 

-55.452 


-1 .728 

17.486 

-43,714 


-2,304 

1 .2855 

- 3.2138 



Contd 
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Table 4,4 (Continued') ; 


x-coordinates 

mm 

Px = 

-0^12 N/ram^ 


V 

1.008 

0.079619 

0.19905 

0.864 

-8.2233 

-20.558 

0.576 

-13.719 

-34.298 

0.288 

-18.222 

-45.555 

-0.288 

-12.717 

-56.984 

-0.576 

- 2.9467 

-59.396 

-0.864 

1 .9699 

-59.038 

-1.152 

6.8797 

-57.534 

-1.44 

18.91 

-52,35 

-1.728 

17.464 

-43.66 

-2,016 

13.535 

-33.837 

-2.3040 

6.5141 

-16.285 

-2.4408 

0.29756 

- 0,74390 


x-coordinates 

Pjj = -0.14 N/ran^ 


t 

t 


X 

y 


0.8928 

-0.35573 

-0.8893' 

0.576 

-12.416 

-31,040 

0.288 

-16.592 

-41 ,480 

-0.288 

-16.964 

-55.35 

-0.576 

- 4.9401 

-59.126 

-0.864 

- 0.17571 

-59.391 

-1,152 

4.9587 

-58.477 

-1 .44 

13.267 

-55.070 

-1 .728 

19.215 

-48.037 

-2.016 

15.331 

-38.327 

-2.304 

11.307 

-28.268 

-2,5848 

-le i 866 

2.9665 
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Table 4,5 : PMMA Cylinder on Steel Foundation 

(Plane Stress) Vertical Load = 125,5 N/mm. Tract- 
ion distributions on the foundation 


X -CO ordinates 
on foundation 
(mm) 

Px = 

-0,01 N/mm^ 



1.512 

0.12969 

0,32423 

1.152 

-13,005 

-32.512 

0,864 

-15,826 

-39,564 

0.516 

-19.343 

-48,357 

0,288 

- 5.1569 

-53,562 

-0.288 

2,8598 

-54,261 

-0.576 

13.075 

-50,948 

-0,864 

17,683 

-44,208 

-1.152 

13.972 

-34.930 

-1,44 

9.5431 

-23.858 

-1.656 

-0,56135 

1.4034 


x-coordinates 
on foundation 
( mm) 

Pjj = - 

5 

0,045 N/mm 



1 .2672 

0,20488 

0.5122 

0.9504 

-11.833 

-29.583 

0,6336 

-15,576 

-38,94 

0.3168 

-16,792 

-49.097 

-0,3168 

- 0,044009 

-55,066 

-0.6336 

5.2803 

-53,566 

-0.9504 

17,779 

-48,637 

-1,2672 

15.733 

-39.332 

-1.584 

11,897 

-29.742 

-1,9008 

-0.01165 

0,027913 


Continued, 
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Table 4>5 (Contimied^ t 


x-coordinates 
on foundation 
( ram) 

Px = 

-0.06 N/mm^ 



1 .152 

0.028251 

0.070628 

0.864 

-11.311 

-28.276 

0.576 

-14,836 

-37.089 

0.288 

-1 8,925 

-47.313 

-0.288 

- 2.6371 

-54,719 

-0.576 

2.7239 

-54,14 

-0,864 

8.9688 

-52.192 

-1.152 

18.528 

-46.321 

-1.144 

14.875 

-37,188 

-1 .728 

11.191 

-27,978 

-2,016 

- 0.4367 

1 .0918 


x-coordinates 
on foundation 
(mm) 

Px = 

-0.075 N/mm^ 


V 

0.9936 

-1,0260 

-2.565 

0.79488 

-9.7082 

-24.271 

0.59616 

-12.853 

-32,132 

0,39744 

-15.575 

-38*937 

0.19872 

-1 8,426 

-46,065 

-0.19872 

- 8.1476 

-53,054 

-0.39744 

- 2,9691 

-54,297 

-0,59616 

0,5159 

-54.535 

-0,79488 

4.2568 

-53,941 

-0.9936 

7.8747 

-52,692 

-1,19232 

16.781 

-49.251 

-1.391 04 

17.449 

-43,623 

-1.58976 

15,047 

-37.617 

-1.78848 

12.108 

-30.269 

-1 .98720 

9.0397 

-22.599 

-2.18592 

-0.75849 

1,8961 


Continued 




0,8784 

- 1.0314 

- 2.5785 

0.576 

- 11,819 

- 29,546 

0.288 

- 15,960 

- 39,899 

- 0,288 

- 12.198 

- 52.017 

- 0.576 

- 2.1875 

- 54,397 

- 0.864 

2.5601 

- 53.938 

- 1.152 

9^,1389 

- 51.665 

- 1 ,44 

18^,198 

- 45,496 

- 1,728 

14,430 

- 36,076 

- 2.016 

10.371 

- 25,928 

- 2,268 

- 0,61621 

1.4505 
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Table 4,6 : Same Material for Cylinder and Foundation 

(Plane Strain), Vertical Load = 32,84 
N/mm, Traction distributions on the 


2 foundation 

-6 N/mm , |jt= 0,01,y = 0,35, E=4000,fiiPa, Radius=31,4 nan 


x-coordinates t 

on foundation ^,,^2 

(mm) 



Applied 
Horizontal 
Traction (P„) 
2 * 
N/nan^ 


0.7536 

0,5652 

0.5024 

0,3768 

0,2512 

0,1256 

-0.1256 

-0,2512 

-0,3768 

-0,5024 

- 0,6280 

-0,7536 


0, 2054331 E-02 -0.2054331 Ef 00 

0.1944955E+00 -0,1944955E{-02 

0.1 90531 3E+00 -0.2004528Ef02 

0,7110624E-01 -0,2364964Ef02 

0,4838224E-01 -0.2575839Ef 02 

-0,1167128E-G2 -0.272291 8 Ei-02 

-0,161 0054E-01 -0,2748651 Ef 02 

-0,6952741 E-01 -0.2628408E+02 

-0.l421907Ef00 -0,2443732Ef02 

-0,2134397Ef00 -0,21 34397E+02 

-0.1732978Ef00 -0,1732978E+02 

-0.532391 9E-01 -0.532391 9Ei- 01 


0,7536 0.45591 13E-01 0.45591 1 3Ef 01 

0.6280 0.1467963Ei-00 -0.1 46796 3E+ 02 

0.5024 0,1450045Ef00 -0,191 31 07Ef02 

0,3768 0.5829449E-01 - 0.2 306485 Ef02 

0.2512 0,4032738E-01 -0, 2546045 Ef02 

0.1256 -0.5558744E-02 -0.271 3859B^-02 -0.01 

-0.1256 -0.21 431 95E-01 -0.2781 172Ei-02 

-0,2512 -0.7800949E-01 -0,26841 36E+02 

-0,3768 -0,1758459Ef00 - 0 . 252681 2Ei-02 

-0.5024 -0,2271789Ei-00 -0,2271789E+02 

-0.6280 -0.1861081E+0G -Q.l 861 081 Ef 02 

-0.7536 -0,1 406931 E+00 -O^ 14069 31 Et- 02 

-0.8792 0.6769729E-G1 0, 67697 29 Ei- 01 


0.5495 

0,5024 

0.3768 

0,2512 

0.1256 

-0,1256 

-0,2512 

-0.3768 

-0.5024 

- 0.6280 

-0.7536 

-0.8792 

-0.10048 


0,2396421 E-01 0.2396421 E^- 01 

-0.4346475E-01 -0.4346475E+01 

-0,1 624831 Bf 00 -0,1 624831 E+ 02 

-G,1627583E+00 -0.2065658E+02 

-0,1 000938Ef 00 -0.241 0768Ei-02 

-0,1573599EfO0 -0.27421 26 E+02 -0.05 

-0,2759253Ef00 -0,.2759253E}-02 

-0,2728'='39E+00 -0,2728539E4-02 

-0.261 4240E+00 -0.2614240Ef-02 

-0.241 3473E1-00 -0.241 3473E+02 

-0.2G85877E+OG -0,2085877Ef02 

-0, 1674826E+00 -0,1 674826Et-02 

-0.312361 9E-01 -0.31 236 19E+01 
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Table 4,7 ; Teflon Cylinder on a Steel Foundation (Plane Stre 

Vertical Load = 75,3 N/mm, Traction Distribution 
20'^ Foundation 

Py = -6 N/mm = 0,4, y = 0,35, = 0,59GPa, Radius = 72 


x-coordinates 
on foundation 
( mm) 


0,1378058E+01 
0.1181193E+01 
0.9843272E+00 
0.787461 8Ef 00 
0,5905963E^-00 
0,3937309E400 
0.1968654E+00 
-0,1968654E+00 
-Q, 39 37309 Ef 00 
-0.,5905963Et-00 
-0.787461 8E+ 00 
-0.9843272Ef00 
-0..1181193E+01 
-0.1574924Et-01 



Applied 
horizontal 
traction (P 

N/mm^ 


-0. 3444666 E+00 
-0.2817621 E+ 01 
-0,371 9509 E+01 
-0.4525787E+01 
-0.5255468E+01 
-0.3239295E1-01 
-0,1268908Ef01 
0,5891 16OE4-OO 
0,1 961 679 Bf 01 
0.4445024E4-01 
0.4974437 E+01 
0,4253855 E+01 
0,3372676 E+01 
-0,7954376E+00 


-0.861 1666E+00 
-0.7044052E+01 
-0.9298772E+01 
-0.1131447EI-02 
-0,131 3867E+02 
-0.1451671 E+02 
-0,1 51 4729 E+ 02 
-0'^,1529499E+02 
-0,1 497321 E+02 
-0.1405972E+02 
-0.1 243609 E+02 
-0.1063464E+02 
-0.8431 690E+01 
0.1988594E+01 


- 0.01 


0.9843272Eh00 
0,787461 8E+ 00 
0.5905963E+00 
0,39 37309 E+00 
0,1 968654 00 
-0,1968654E+00 
-0,39 37309 E+00 
-0.5905963E+00 
-0,787461 8Ei- 00 
-0.9843272E+00 
-0,1181193E+01 
-0.1378058E+01 
-0,1 574924 E+01 
-0,1 771 789 E+01 
-0.1968654E+01 


-0,1 365928 E+01 
-0,3160373E+01 
-0,403031 3E+01 
-0,4747751 E+01 
-0,5016703E+01 
-0.1055871 E+01 
0,61 371 92E-02 
0.1116238E+01 
0,2189050E+01 
0.50441 95 E+01 
0,47651 30E+01 
0,4029027E+01 
0,3078984 E+01 
0,2094911 E+01 
-0,1726547E+01 


-0.341 4821 El- 01 
-0.7900933E+01 
-0.1007578E+02 
-0,1189438E+02 
-Oil 385351 E+02 
-0,1 535221 E+02 
-0,1534935E+02 
-0,1519338E+02 
- 0.1 482261 E+02 
-0.1369607E+02 
-0.1191282E+02 
-0.1007257E+02 
-0.7697459&-01 
-0.5237277 El-01 
0,431 6368E+01 


-0.07 


0,787461 8E+00 
0,5905963E+0G 
0,3937309E+00 
0,1968654B+00 
-0.1968654E+00 
-0. 3937309 E+OO 
-0.5905963E+00 
-0.787461 8E+00 
-0o9843272E+00 
-0.1181193E+01 
-0.1378058E+01 
-0„ 1 574924E+01 
-0,1771789E+01 
-0.1968654E+01 
.n, 21921 34E+Q1 


-O ,1 666333E+01 
-0,3274758E+01 
-0.4111599E+01 
-0,499701 0E+01 
-0.3113656E+01 
-0,1080090E+01 
-0.2692308E-01 
0.1 117485 E+01 
0,21 89340 E+01 
0.4945521 E+01 
0.4795838E+01 
0.4061 970E+01 
0.3130260E+01 
0,221 3079 E+01 
-n 1.5Q430E+01 


-0,4165833E+01 

-0.8186895E+01 

-0.1027900E+02 

-0.1249253E+02 

-0.1479780E+02 

-0.1 536889 E+02 

-0.1542662E+02 

-0.1525008E+02 

-0.1486604E+02 

-0.1375445E+02 

-0,1198960E+02 

-0.1015493E+02 

-0,7825649E+01 

-0,5532697E+01 


-0.1 
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CHAPTER 5 
CONCLUSIONS 

Based on the results and discussions presented in 

Chapter 4, following conclusions are made: 

(i) In BEM displacements and tractions are the variables 
whereas in FEM displacements and forces are the 
variables and this fact makes BEM more suitable for 
contact problems, 

( ii) Linear eluents give better results than constant 

elements in case of traction fields as well as displace- 
ment fields. Moreover if the load is applied in 
increments and the discretization is done with linear 
elements, introducii^ previous step's displacement fields 
in the existing model is easier, 

( iii) It is observed that solving the contact problem in 
natural coordinates system is preferable as it can take 
care of any configuration of the contact region easily. 
When large constact region is present, tedious matrix 
operations are needed if cartesian coordinates are 
used, but in a natural coordinate system this problem 
will not arise. 

(iv) While all primary solutions, displacement fields and 
traction fields, bdiave well the secondary solution 
(moment about the centre in rolling contact) show 


instability. 
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Scope for further work; 


(i) To get the secondary solution, moment, in an accurate 

manner alternate ways like higher precision calculation 
and introducing moment as a degree of freedom may be 
explored. 


(ii) To simulate rolling wheel at all speeds boundary 

conditions based on ’creep* factor (if two vdieels are 

in contact and their velocities in the contact region 

2(V. - V^) 

are and ¥2 then creep, | , = ' y ‘ 

be used. 
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appendix a 

Boundary stresses can be calculated by a method 
presented by Rizzo and Shippy [17] . This method 
obviates the necessity to evaluate singular integrals. 

The method is described below. 

There are seven unknowns on the boundary, viz, 3 
stress components and 4 displacement gradient components 

seven equations are available all these unknowns 
can be evaluated, o. = cr.. n. ( n 

‘ j 1 j 1 ‘ ^ 

supplies two equations #iile stress-strain relationships 
supply three equations. Two more, equations can be obtained 
by evaluating duj^/ds numerically along the boundary and 
equating it to analytical expression for the same 

du. dUj_ dx, du. dX 2 

■ar = di^ Hr + 

= “1,2 "l “ “i,1 "2 

where 

n^ , 1^2 - normal components in n.j and n 2 directions 
s - indicates boundary. 

Along with Eqn, 2 there are seven equations for seven 
unknowns. Solving these linear algebraic equations stresses 
on the boundary can be obtained. 
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1 ( ms tari function oolnts ; 
a', joimt ft the corresponding pts. 

K mr I'lsssian pta. and twelve ' gaussian points are given 
(,) ...'Vd'r ddly id qaussian oolnts are used in the wori< 


‘.•ii’0CU=y.yhii3h) ; 
iPPCD^-o.Bhiuei 
m i )=(). 34/8548 ? 

wTt4)=0, 3478548 ? 


GFP(2)=0. 3399810 
? GFPU)=-0. 3399810 

6»T(2)=0, 6521452 
WTCSlsO. 6521452 


I )=t’. 1252334085 
Uf"-'12)=-GFP(1 J 
Gi v-l i)-( , 10/8314989982 

CP j) 

gr-'M »)=n. 587 3179543 
1"1) = «gfP(.5J 
GfiH / /=' .7h99u2b /42 
t.,r,>(u Js-GFPf/J 
Gi'ot'Msu, 90411726637 

Gi'of. i j )=u,9)a560d342467 


MJ 





ViT<3)s 
WT(4)= 
?WT(53= 
;wl'C0 = 

,*WT(8) = 
,‘WTf9) = 

jwrcu) 


0,2491470458 

wl’cn 

0.23349253654 
WTC 3 ) 

0.203167426723 

WT{5) 

0.160078328543 

WT(7) 

0.106939325995 
sWT C 9 ) 

=0.04/1753363865 

;^X(12)=WTtll) 


P iii'ji;; 'FjrifcRliiu GHhA’i’X' , , . . 

^vdtein matrices IHJ & [G3 are calculated 
c. .’.MhATXi Xu,yn,) 

III * n;j are stored j.n a tile for reusing as long as the 
r- wr Uiritos of the nodes do not change. If coordintes change 
i v'u *rusii IHi [Gj will he written in the same file 


m 2/ 
. ( 2 I ) 
h ( /; / ) 




'p.\usr: 'E6Tl?!UNG ^ rwTAfT 

nlinost all the wor<< is fS^pfaki 

,.4 j F.iihrTC’MXANGrKCOdTfXCJ ,yN,Xl, iI»C0 »ll/E0 ,PEAll/ 


GT* t.uO 
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'■'uLo rojrine introduces the contact conditions ana rearranges 
t'to '-itrlces in sucn a way that finally UJ iXJ s [BJ form is 
ot'txinnri an 1 tnls is solved by LU cJecomDOsition method. This 
rnuti'n; »lso ifi.Jintains tne iterations till needed. 

HU.n.n I au;: Cf^lACT(TOTANG,KCONT,XN,XN,XX/n,CONLEW,PEAK) 

u rVu i jcn J A'aiupl) ,XU( JLNPl ) ,KC0NT(2,32) ,M0(2) .C0f^TATC2rid) 

Oa’II" Idl.) i K0bdLPt2,20),l’Jl)SLlP(2,16),CnHP(2,lb) 

UI-i.:'j;Vi!M T0X4 ■•Gt2, Jt>) ,DEGl C2),DEG2(2) ,X1(20) ,XT(20) 

DI A'-.urui} 1 urjCA'JlC5U),GUCA12(503,Rn(2»l6} 

CD iO- »■.;/ i,\T/ dM C 1 9i , 19 3 ) , GM C 1 9i , 1 93) ,MiSV (193,1) 

C'i - 4 i.J/£P>ii;/GFP(i2J ,WT(32) 

Ct) i.uiu/f:u3f.lT/G£l ,POn ;CM0b(2) ,PRC2) ,AMU,Pi 

CM 't .J’J/hCS/QldHitl93) , lHACT(l93),BV(l93),KUDE(19i) 

Cu "i.’iUU/PPGi/Ab/ AtlM,N0C0K,lPm)B,tlC(2) ,NbU8, 

1 l^jUH(io),iNDE(lO),COWTlr) 

Cinh.i J/T<<AhS/.'iuO'lPt2) ,9QUnViU(2) ,AL«)CAT(2) ,}lANGt2) ,SRARCH(2) 
Ct3'ViJf^/i)t I/IU, I.H2, IHM, 1HM2, INPl 
Cri't.'l(Jh/SUHCri/P(3 9i) ,WKSC193) 


M.r./ 

(iC=G'i.)hC t ) 

NHCFlX = 2#ii''i 
hBCSHdCFIX 

the followiny Input helps starting the program from an 
intermediate stage , 

READt23 ,*)lTEh,UU( 1 ) ,NU(2) ,i[iODLEF,NODRIT,lRFIX 


DO 820 Isi.f^OCON 

CaNTAT{l,l)=TOTAhG( 1 r I > -tOTANGC I , I”! ) 
COMTAT(2jl)=TUTAhG(2,I)-TOTANG(2rI*-l) 
CONtlNDE ■ 


Beginning of iterative loop 


Kiuf:=<) 

ivT , : ''5 2? 

(2/) H I 
''M.-'i-r/j g-4 


i M ( i ) : i 0 


. of noces in contact on the left side of 
axis ot symmetry and ho{ 2J similarly on the 


?-,i, j,(t side 


• ),duC3) 


i. : -'’.li is total contact length. 


: = , l.i‘:corit:l,2*»OtU-U)) + A!3S(Xu(KCUNTC2r2=!‘NDt2)-13)3 

jI I icr4u'^^:es tne site of tratrlces LHlU & [GMl from 2R x 2H 
t) (2. + 3) '< (20 + 3) 

Mi ^ i, jCoLuCATI ,LOCAi2,KCDM'fXl'i,Yti,NTOT,R,NO) 
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f-'”! ’.i jf 1 ~ I If 'i't'O'X 


<’ j-.iiirr' [X 

it* 

if't iHPlA.iiO.*- j!l!j(2J='} 

-fh..' »;oi.lJvp(j loou C 6^ DO ... CONTINUE )introduces 
tnc contact coivditions . 

itn 6*1 'iU.1 = i,lT£R 

Ur 66 K=!l,2 

ir (6.60.1 .AwD. Li* IX. EQ.O)GO TO 68 
LF(.'v.Cy.i^lO(j lu 78 
it INfM.GT.fiUtl ) jOU 10 6S 

CU I'p 6 1 

t'U I 1= i;ja 1 + 1 

(i U ‘I f ) 61 


IPCTKFI X,Ed,6 )c,0 TO d9 
IF( ^U'i.01.iv0t2))G0 10 65 
oO Ti.f f>i 

ftQ(2J3:yu(21 + l 

cont iiiiji: 


1T = I'!IIM-1 

FHICsAHU 

NDSUP nas toe nun.Ders ot toe nodes tnat are slipping 


iFCliUSOIPCK.NUM) .eo.KCUNTCK,2»CNUM-l )+2)lGD 
//SD1S=(XN(KC00T(K,2*NUM-1 ) )+CONLEn 5+ZETA// 
//iFtITRADD.EO.O) SDIS=0.// 


TO 575 


DO 650 I = l,N'rOT 

G.ICI ,2*KCUftT(K,2*NUMl-l)s: 

1 Gni,2*KCUNt(K,2t^NrtM)-l )-GM 1 1 , 2^KCaNT (K , 2 ’•= ( NUM-l ) + 1 3 -1 ) 

i1CX,2*KCUNT(K,2’i^M!M)-n = 

1 I i( I,2*KCUNT(K,2*Nl:w)-l) + HVI(I ,2*KCDNT(K,2=»'(MUM-i) + U-l) 

U (' 111! 1 -101,102 

> n i , i 1 = ' e V 1 1 , 1 } + Hrt { 1 , 2* ftCuMT (K , 2 * ( N UM-1 3 + 13-1) »DEliSl.P 

. 'I ■ , 1 ■ 4 

•I I, : ,'t c'Ui,13 + H?UI,2 + KCOi>.-l (K,2 + {NUM-13+13-13*DELSLP 


1 

1 


C 


Hi , If+CO il (n,2 + ru,j6) 3 = 

. ( i , 1 * **. C 0 *, Ti K , 2 + N u 3 3 -Gi<i ( i , 2 ♦ KC ON 1 C K , 2 * t tiu M- 1 3 + 1 3 3 
■;*U ,(.,})= »HytX,13 - HMf l ,2«=KCOfiT(K,2»tNOM-U+l)3* 

(. f ; . ( ru, ; u i T ( !*> , 2 * ( f , U N ■» 1 3 + 2 3 ) - Y in K C 0 N T ( K , 2 + ( N U M - 1 3 + 1 3 3 3 

1 (!,.. + K C u t T O , 2 4 Mi ?- 3 3 + N +i C 1 , 2 * KC: 0 N T C K , 2 * ( N U M - i 3 + 1 3 3 

.‘IT-.'L. 


.i I ■ ! : . . 1 1 (. i . .1 . r. + ; 3 = 2 ^ K.c a . t ( k , 2 » c n o +> - j. 3 + 1 ) - 1 
■il i. ■, > 6t-6C ( 1 J 3 = /*r;CC. ^T (K ,?*NUr 3-1 

. i i U . . 'I : ( K I . C + 23 = ♦ R c t‘ <’ ■’ T ( K , 2 * C & I N'1 - 1 3 + 1 3 
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3 

i 

C{ 


: - [ ‘ : r, ' 1 3; U‘ X . : C + 2 ) = 2 » K C 0 N T t f! , J u J 

!>u f>:yj i=i f fi'j'OT 

t iCl CRUK-n + D-l J: 

- •■»’< (. I f 2 * KCCiS;f t K , 2 ♦ KUf? ) - 1 ) »GE 
H •■ ) C i , 2 * K, (• 0 N ' r C K , 2 ♦ 1 U F- n + 1 } J = 

Cn'iXl’;tiE 

. J ' « C ~ i “i C + 2 
•iVj bb‘~j 


,M,/iNI).fdfF,l':U.NCtK)3GQ TCi 411 
,lr A'’D,NUM,i;y,NO(K) )GO TO 411 

GO T'-i iy'/ 

COmTXI'MJE 

t iaC=COHf'(K,f)QtK)~l )*AMU 
V4;1'U. (.b , *J ^-yXC 
Gij Tu iy(j 
CGb'Uiiiji: 

X r ( T f.C T 1 2 ♦ KC U 0 1 1 K , 2 » ( N y H - 1) + 2 ) - 1 ) >^TB AC T 1 2 ♦ KCO NT 
(f' , 2*^{rH)M-l )+2)) .GT.O, JFRIC=-AMU 

WP:ITt,(5 ,F) F hlC 

on yyi- i = i,ntoi 

iF(')U-l,Ey.NO(K3 jyK=Na(K)-l 
G UI#24KCUNT(K,2*NUM)-1)=: 

1 G'UI,2»i'vCUf*T(K,2»NUM)-lJ+GMtI,2*KCONTCK,2*NUM3)/ 

2 («*AFtU*CUfiPCK,JK)) 

GM(1 ,24KCUyT(K,24NUM)-13= 

1 GM(I,2*KCUf4T(K,24NUH}-.U-GMH,2»KCONTtK,2*CNUM-13 + l)-l) 

GM(I,2>»KeUNT(N,2»NUM)-l)= ' 

1 GH(i,2FKCONTtK,2*NUF3-l) + G«Cl,24KCUNX{K,2*tNUH-U + l))/ 

2 (AMU^CGMPCK/JK)) • . 

MBVtl,l)= HBVCiJ) - HMCI,2*KCnMTCK,2*CNyM-i) + in* 

1 (XN(KCUNTCK,2HNUM«l) + 2)3-1£rKKCDNT(KV2*(NUM-l) + l))) 

HM1I.2»KCUNT(K,2*RUM))= 

1 HHCL 2 FrXONTXK, 2 »NUM)) + HMCl, 2 'i'hCrjN'i‘CK, 2 »tHUM-l) 4 l)) 

HMCI,2»KCONTtK,2*CNyM-l)-»-U)= -CM C X , 2^KCaNT C K , 2»NU(4 ) - U * 

1 GE ' ■ 

CONTINUE 




CAT! (KXNC+1)=2»KC0NTCK,2»(M0M-1 )+l ) 


y n c A T 2 (, K I N C + 1 ) = 2 * K C 0 N T C K , 2 * U u M ) - 1 


1 


C() 

: i, ■ ( 0 , 1^3 (. 0 .'C-VfUl),I=:t,EXWC) 

: i.'. '.r.,-- )U', .C-.,!;itX3 ,T = 1 ,Fi;.c) 

A. ■'■‘r If xc J '■.;‘ir;;).3ry conditions are lOfjlimented oelow 
5.' ( 1 A 
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tiiue of lAl iX^ - i B > is found out 

c<,bi, 

= ii'i’nr 

Cabi. MATIMC'i.'i,b8C3,'^BV,l,Dfc;T) 

C -Abb t'r,4APt*' t Hr>) , I UM , 'iDV, NBC3 ,SOL, WKS, IFAIL) 
i)r, 2* d 1 = 1,UBCJ 
MuwCi ,n = suj,t 1,1 J 

Cb - ; '.ut/i., 

Souuciuij ot (AJ i A '«■ = i B 1 By BU decomposition 


CALL F <» 1 , i T F C li QC J # H B , UPI , F , D P , I F A 1 L ) 

TYFB ♦,lFAiL,DF 

C7»LL Fd4Ay F C. NHC3 r 1 , HM , IHM , P , HBV , IHM , IFAIL ) 
TYPE ♦,IFAiL 


Output below , using variables LOCATl , BDCAT2 KODE values in 
MBV I <,X> ) are properly related to tne corresponding nodes. 


UU 310 K=1,2*NN 

lid J=:l ,KINC 

lMi.,.:i,L<iCAH(J))GO TO f,?0 

’(b.u u ‘->,125,130 

■ rv.>I.U' )=b<# IK) 

■S (.K» 1. )1'GE 

^r' I.J 

■ 1 ;j . ’ )(.■.,*“ 1 . . V (. n , i ) 

1 . '■"!•. .1“-; 

I 1 f'i 


t .bU. LtlCAT21d)/2*2> GU Tfl f,71 

'.r;\C'UuCCAT2(J)) = MBVtLUCAtUj),!) » GE 
Lf; lu in 2 

TiLil-UiiUCAt2td) ) = BBV (LOCATl * GE 
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(.! u ’ * i' I n u £ 


Cri nfl nuti 


tne tractlofis and displacements of nodes other than 
t;io5'.' t.ouni oat from direct solution witn the helo of contact 
Coalitions. 

n () 7 5 •<> = 1 , 2 
0‘i 67h J = l,NrKK) 

2/J5r)lS=tXhUCur4‘tK;2=«'riUM-i))+CONLEh)*ZETA// 
iF(d.c;y,HutKi)dK=Na(K)-i 


TH act ( 2 * KCUMX ( K , 2 ( J- 1 ) + 1 ) -1 J = -I'F ACT ( 2 ♦KCHN'r ( K , 2 ♦J ) - 1 ) 
TflACTt2»KCUHT(K,2*tJ-l) + i) )= -TRACI (2^KCUN'f C ti. 2+a ) ) 
|)TF.PH2»KCUNT{K,2*CJ~U + 1) Js DISPLC2=^KCOi'iT(K ,2*J) ) + 
l’!’!lCKCUNT(h,2^Cj-l) + 2) )-rNCKCfJNT(K,2*CJ-i) + l J 3 J 
iFCNf)3lii.PlK,J3 .EQ,KC:0NX{K,2*J) 3 GO TO 680 
bf) '!'!) 685 

TRA<h(.2*KCUHT(K,2=»J)3 = CTRACTC2*KCONTtK,2<^JJ-n/ 

<-AMU*C0MPtK,JR)3) 

l'RAC‘l‘r2*KCUNT(K,2’»CJ-13 + l)3= -TRACT C2*KCOMTCR,2’^J) ) 

GO TO 675 

0ISPL(;2^'KC:UilT(K,2^tJ-l 3 + 13-l3=OISPli(2’t'KCUNTCK,2*J)-l) 
1F1NM31C4,164,105 

i)ISFo(2»KCONTtK,2*(J-l) + 13-l) = 

Oi6PL(2»KCaNT(K,2»ld-13 + 13-l3-OEXiSliP 
GO to soej 

i)tSPU(.2*KCUNT(K.2»(d-l) + l 3-13 = 

DISPL(2»KCONT(K,2»tJ-13+l)-l3-OELSLP 

COHTINUb 


CONTINUE 

Secondary variables viz. Force-x , Moment due to Force-x# , 
Moment due to Force-y are also directly obtained by introducing 
appropriate equations in the matrix system. 

Alternate way is to calculate separately from the results. 


AMX 

FI 

AMI 


HBVi NBC 3-2,1 
M«V(1«»BC3-1,1 
MBVtNBCirlJ 



i (24,7073 
i- OK .:,i (/,lox, 

1 1"X, 

2 li'\, 


amx.ami.fi 

■'moment due to x-fqrce = 

•moment ,, ,, T-FORCE = 
'Y-ffJKCE IN CONTACT AREA = 


f 



f 


E2(},8,//, 

E20,8,//, 

E20.8,///3 


I oil' .ri.i I Stores important input & output values 

tiji fietare ise in otner programs. 

>u;:-iT ■ c Pt i 

-i ■ I ' ' n * 2 1 ; ) /I > « , i N , K C , M I , ( 1 C H , X 1 , H , T H A C T , D 1 s p h 

it' ( ''.nu.lj nO TO JOB 

c; a 

C:V’p .stores tne Kl-.jti ot AMU in each slipping node, 

Cr iter in is Tx + a iu *-'iy = ii 


UO .joii Ksi,2 

::»>> iUa j=il ,HU(G) 
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ir'( 1 

U(! 


AC’CCi'-f'KCi/fjlHK 
Cf) liniv,.J J = 1 
TO J9h 


I J-l)»TRACT{2!f'KCUNT(K,l) J J.GT.y. JGO 


fO 399 


U.'! '•'!/ iK'i 
UP .1' 7 U=l,2 

Ut' 107 J=lrh(jtKJ 

IKUUty.l JCflPP{K,d3= 1. 
iF(K,t:y.23Cf.)f:P(K,J3 = -l, 

Crifprif.Mih 

J'O C t l,,l)15PLC2*I-l) ,Dl3PLC2*I),TRACT(2*i-13 ,TRACTt2*I) 
1 J , 3 S'P , 4f.» ) 

'"i I ‘ X L, (5 7 7/) 

(••iij-in/V'J'Ux, ' 'ilfPE IN ZE’HQ value ' ) 

Accept », .u.ip 

xrtiuuxu.Ed.nuu tu 544 

irCTHACTCi’^KCuNTa fZ^CNQtU-U+l )).GT.O, 

1 .na.AiuSt r!U\CT(2*KCUMT(.l ) + i) )).Ll'.AiilM)GU I'O 541 

if LFIX = 0 lett side contact region (w,r,t vertical 
axis of syr-iTietrv} can grow . It LFix = 1 no more growth 
ailowQi . Gimiiariy tor right side IRFXX is used for 
controlling growth. 

Ih CIPHTUjU J2,y.l2,93i 
iiFlXsp 
GO to 542 

fiLlCl )=UO(l J+1 

noolef=uodlef+i 

GO TO 542 
0FIX=1 t IFHIL=:1 

I F ( TN AC T ( 2 =^KCUNT C 2 r 2* ( NO { 2 1 - 1 3 4 1 ) 5 , GT . 0 , 

1 .fia,r.BSCTRACTt2*KCONT{2^2=f=lfjUC2)-13 + l))).IiT,AHM3GO TO 543 

Lf'lTPHIHjqiO, 930^931 

Gi TO 5<H 

c. r2} = L(iC2) + l 

> 1 I'Pils'tOLipl'i. + i 

I ■ 5 i‘l ' 


1 

1 




' .'..•L. -.,1. 

■ J = X ' r A . 

■"l/.'Xuiia 

3 A f X r ^ X ^ 


' , 1- X , • U I GPLACEhELT X ' 
'I’MCTin.. V',//) 
r u 1 A P H 2 4 1 - 1 ) , D 1 5 P L c 2 

hl2.5,7X,Ei2.5,7X,fc.l2 


/5X, ‘UlSPLACEhENT Y',5X, 

♦I) ,TKACT(2»T-1) ,TRACT(2#X3 
.5,3X,E12.53 



10 




CO''JTi/jUE 
•,'U li% K = 

iJU JJo Js:i,Nn(K.) 

'(T (. K , 2*j > -1 3 -DiSPL C 2 + KCONT £ K , 2 ’»' ( J-U + 1 ) -I ) 
n /) U^fKK,aKE),JKL = J ,nO(KJ) 

H/, i.".). *(:;;'iTi.%«iNG slip* 

C.'il.J. r'iilPCKCU.'^T.MOSEIp^KUDSLP, lADD, 

_ ^ .JAJ)0,KADD, lTHAUU,riEH,KO,D&LSLP,MSLIP,IKi , iK21 

I'nr ^itatlc contact proDleirs activate tne following 2 statements 

•inis I p =5 0 , 

■^o to 20J 


<‘'.5 


iijf A 


r 


leips In controlling tne a subioop ot iteration to 


introsuce tne relative displacement' in no-slip region also 
to s if! J late rolling contact conditions. 


"■-ni,r.p=i 

lFUi.ti.fc-0.0 .ANU, aK2.GT,0))MSLiP=0 
.‘.f uiM.Fi|<!i.n-DnL.U.P 

<l:U'ht^ril2 ) ae.LOlt .RRDIF 

Fnl' <<ATCnx, 'REtniF s % K 16 .5 #/ , 5X , * RHDIF = ' ,r,16,5) 

IE’ lABa(.HKlUEJ.GF,U,OOOal)GO TU 5 

■13MF- ■’ r HF:L=n, ; ph=o. 

CAii!. 6LXRCAC0 nunn6ljIP,KQDSLP,IADD, 

, J A n n , << aD j , I T RA 00 , I iEB , NO , DEL5L P , MS n I F , 1 K 1 , IK 2 ) 

GO TO 20J 

continue: 

REOsUtiLfiOP 

HRSREOUIE 

WHITEC5r*K C(N06 0lplK,U,KGDSLPtK,I) ) ,1 = 1 ,1TFR),K=;1,2J 

FORMATC/,5X,bE14.4J 
TIPE lADD 
TtPE. *, MSUP 

PAOSO 'AFTEIH ASSESSING SLIP' 
ll’(IAnu.f:g.O)GQ TO 860 


i.j i "I 


llib 


f/. , 


• a , i i J “t j ; iLL 2 =KCOPT 1 1 , 2 * (no t u-i )-i 
‘■J.’i.'-lU lP.h2 = KCOnTC2 f2=ECN0t2)-l)-l 
I .■» ! ’ fr. ; ,T! AC. f f 24 TLF ) 'T^hr > .TRACT 


i. '.I' , 

.0, 4 j 

•-= i 

= ' -O.-i j 


. . , , . M c ; I C. 2 4 1 6 H ) , T I ? A C T ( 2 ^ 1 R E ) , T R A C T 1 2 I L E 2 ) 

' >\r 1 1 . :.'4 iLi :;- 1 ) , tract c ? * l HE-1 ) 


Oii .i u 1 1»'.; 2 , U'i .J ) , K 


r m T ii A C ■ S: C 2 * I L E > . G T . T H A C T C 2 ^ X 6 E 2 ) ) N U ( 1 ) = N 0 ( I ) - 1 



it 


>iO TU bgl 

1 1 .’ C 'j : ] < A c T (, 2 » 1 H E J , G T , T H A C T ( 2 ♦ I K R 2 } > N 0 C 2 ) = N □ ( 2 ) - 1 

'I'.'UlUi-, 


t-(:Af. 3!A;F.Ml:,.^4U(l J.rH<,N0l)«tT.NB,NDi2))G0 Tu 155 

I 1 1 I. !<AC; T ( 2 » X i, L J } , Gfc . G . . AND , TP AC T C 2 ♦ X RE ) . GE , 0 , J GU TO 153 
■’Ut.iGlTHnCTU^'IUEJ J,EE.ALlM,AHD.THACTC2*iHSj.GE.O.)GO TO 153 
-i<AC'l(2»T.GL)KGE,G..AND,A.BS(THACTC2l'XR£)J.LE.AGIM)GO TO 153 
L- t,,!'-GU’SGiCT(.2"+XUlJ J .LE.AGIM.AND, 

1 AH5iTHACT(2*lHh}) .LE,ALIM)GO TO 153 

IFUUUlC.Ey.l )GQ TO 153 

■ l'i cF=l.TEHHt I 

TYp!: iter 

IF tX't OU.lFX.'-UCUfj+l J GO TO 3RB 

J.i’CGF LX »EO.O>ISirJCU=:NO(l )4l 

] F(lHFLX,Ey,i>)r,0C2)=fi0(2)+l 
TYFF4, LFIX,iPFlX 
GO T1i 155 


XFLALU>LTHACTl2^XXiEl).LT.ALlM.AND. 

■AboCTf:ACTC2*XREj5 .LT.ALiMlGlJ TO 540 
iFtXTRAOii.LQ.OjGO TO 540 

CAI 4 L /.OOUCCAliXn ,,iUL2) f T0TANG,C0NTAT,IND,INDIC»X«, YN,R,ITRADD, 
•DEGl ,OEG2,K,C0RTX 

CAIjL G(iHATX(XG,Yfl.) 

GO TU 155 
COWTXHUE 

NHsl 

Dn 7^0 K=lr2 

on 700 Jsl.NOCKJ 

lFtK.EQ.UCCJHPLK,0)= 1. 

IFCK.E0.2)CUMPtK,0)s 1, 

Ci.’-. '.ViG .)!, - - 


.CT!2»K(.:LMTtl,2 7) 

t tj Cl-i’jOEN 

^ *!Gi 'iiCl GEOGTH = *,E20.0 ,//) 


:n, ' ,‘> 4 ., 'GXSPLaCEMEOT X ' »5X , 'OISPbACEilEMT Y',5X, 

A* ■''i'‘?A(;TlO:4 Y*,//3 

\ ' ’ 0 ) L L I , j i 6 F L 1 2 » 1 ' U > U 11^ P JU ( 2 ■+ i ) , 1 H A C T X 2 ♦ 1 - 1 ) , T R AC T ( 2 ’t' 1 3 
,,AG,BX,in.:.5 ///,E 3 2.5,7X,E12.5,iX,El2,5 3 
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f* 'li .i ,rt IMX# 'CijtjTACT ELtMENl'g NOT ENOUGH ???????') 

Hi r,'<t 



m mm - m • m mm ^ m * m - mm mm mm mm 




this pro’Jram jets the Input values from CONTACT program are ♦ 

useo in tiis program to calculate stresses at internal points 
as >«reii as t'lstresses at the bounoiary. 


MA[ti program . some of the values have to be hand supplied. 
iOPhiCit REAlf»BCA-ri,Q-2) 

DI -iChSION XF(9 7),YE(97J,X1C20) ,¥1(20 J ,TRAGT(1935 ,DISPLC193) , 
NCt2J,GFPC32J,^TC323 ,IfiDB(10),lNDEa0l 

CO.'l.HUIj/A/GFP,Wt,FI 

CO'-<MU{'i/B/GE,PUl 

P 1=1.3415927 

UPK!: I ihaT=2l,,UEVlCE='DdK',FlliE=' STRESS. IN’) 


NC - ho, ot rujies 

Gc - .Shear modulus of the roller CGE of foundation shouid 
shouict jje Given if stresses in tne foundation have to 
be found out, similarly for POI also) 

POX - Poisson's ratio Qt the roller . 

PEAK - j’eoK value of tractlon-y j used to calculate stresses 
andlvtically when loading conditions are symmetrical, 
tsnin, uy Smith & Liu are used) 

COHLOij - naif contact length ? supplied to EXACT subroutine, 
NUMB - Mo. Of oouoie nodes, 

INDB ” Doubie node beginning 
INDE - Uoubie node ending 


C0WLEN«4, ; PEAK = 9i. i 1NE=96 ? GE=1483.48 / P01=0,35 


OPEN(UNITs28,1)EVICE=*DSK'3 
REWXWD 28 

HEAD 1 28 ) XE , YE, NC , N I , MCH , XI ,¥ 1 , TRACT ,DISPL 


Vil li’? of NDUB, INOB, I. INDE are read from STRESS. IN data file. 

K! ■ i,*).'!DUB 

Ki 1, aMDRCi),INDECl) ,1 = 1 ,NDUB) 


I*' .e 'Loo to ui'/e internai points other than those from 
d:iie, foiin statments have to be used. Otherwise 

t'lf "’ lO.jrtlvated . 


mmmm mm 



^ m m ^ mm ^ m ' m- _ - mm 
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!>ri i. / l.sl ‘.iK 

i i,> ) I »>-'I*>P!K2=i=i-l ),nlSPLC2*l) ,TRACTC2=^I-1) ,THACTC2*1) 

r ».i ; I »,V 3 .i:, 

,'*■£- calculate stresses at internal points 
A.K-M’ linear eiowonts are use^i for discretization. 


*lf 

IP''**" 


IP 


CM,!. 5 !■ -ASL . i I U *1 B , I fJ DB , 1 N UE , X E , YE , X I /X 1 , N I , TR AC T , D I SPL , NC ) 
l iTi IS for internal dlsolacements I 
IM’U H, Is tor internal stresses. \ 


For constant elements, 


CAub iUT!lCai,nE,iiC.XE,SfE,Xl,YI.POIfGe»TBACTiOISPE,MCH) 

CAnu .'.'Jl’Hinuui Ji,n4E/YE,iHAdT,DisPE,!«C,MCH,NE) 

P,)r tvr»e ein.nents analytical results can be obtained 

u$j(i i r,XACi' s inroutine. (limitations are explained In' thesis booR) 

Cam. iMACT(C!.MbH:A,PRAKfXl,^I,NI) 

STnPr’iibn 

Bant. iriA i.E.diHb [s used to calculate the contributions of 
eacb olei.uni In the stress calculation, (For Linear elements) 

S.JiCM.tUTCh':, ELSTRSCNl,N2,XN,YW,XEl,YEljrXE2#YE2. 

1 Dili ,i)2ll ,!iil2,h2i2,!)i22,D222,sfll»S211,Sll2f 5212,5122,3222, 

2 TibAcapuibPb) • 

IMPLICIT BEAb'^Hl A-h,0-Z) 

DlAEUdlO.I Xw(i2).V'^C32),AN0RMC2),FIC2} 

PIMFBhSIiJ.l SFP( 323 ,aTC325 

DlMEH.SIfJiM i)T,.SPu(lPi) ,THACTtiy33 ,BV(193) ,KUDEtl93) 

COdHUB/A/GFBf 

coMHun/ij/Ge,Pui 

ELENG'fssyHTt tXE2-XEl)*»2 + (3(E2-YEl)*»2 3 
ANORM(l3 = (5fE2“'iEl 3/ELEhGT 
ANOHM ( 2 3 B C XE 1 - XE 2 3 / ELENGT 

D1 U=CJ. ;D21 1=0, ; Ull 2=0. ;U2 12=0, ;D122=U.; 0222=0. 

SI 1 1=0, ? 521 1=0. y 51 1 2=0, ?S21 2=0. ?S122=0,; 3222=0. 

CALL GAlj55(XN,YN,XEl,YEl,XE2,YE2rNGPrGFP,WT3 
on 20 K=l,NGP 

F1(13=-0.5»(GFPCK3-1.) 

FT.C2)= 0.5»(GFPtK) + l.) 


,-:=h' i c i j tTHACT(2».Ml- 
! ( 1 J ?T!iArTr'j*i 


_ , , _ +FlC2)=t'TPACrC2*h2-l3 

n'Ul 3 »TraCTC 2?N13+F1 (2)*TPACTC2’I'M2J 
. = ' £ i i 3=*ni>5PiiC2*Nl-l 3+FT (2)*DI5PLf2*N2-l 3 
i V =! ■ T 1 1 3 * D X P L ( 2 1 3 + F I ( 2 3 ISPL C 2 3 


X : I • . - I ( . 2 F X t,, 1 ) / ,' . + GF P ( K 3 * I X E 2 " XE 1 3 / 2 , 0 

,tGFPti<)»CY£2-YEl)/2,0 


(A) 2t (Vh^YVstK) )»=<'23 

(5 


r r; 


k' , J A i I, ' 

r\ j = i ■ i i' .) '» A ) / i '! ’ 
r j-Y .)/,:■ 

.‘i M; 1 i'A )t. 1 Mb.2^-AHUF<h(2) 

E * 'b . / (. b ) 1 b 2 * •> , * B I * ( 1 , -PU I 3 3 
C- l ./(..»'> l -F 1 1 . -BO i, ) -FROIS) 

E~ ' , ""A * £ ■ 

Li . ; = ! ' t l +c ^ t. E x m 1 + 2 , gl » ♦ 3 3 ♦ ( K ) »ELEhGT /2 , ^'Tx 
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1 K f t t,“E J »ik 2 + 2 .*Hl^ 4 ! 2 *R 23 *WTCK)^=ELEMGT/ 2 .*Ty 
12 fC’«th♦H 2 +/,*Rl»=^ 2 *R 2 )*WTCK 3 »ELH^^IGT/ 2 .*TX 
■<r/. ( j/i2+C*(.E^R.l+2.=«‘R2**2*Rl)*WTf K3»ELENGT/'2.»Ty 
uur:= .'j :.::>fL*C(.“f:.)*Ki + 2.*R2»*2*Rl)»WTCK)’!'ELENG'r/2.*Tx 
!»? 1 - -- V'^ 24 -C«(;E»P 2 + 2 ,»R 2 =«'* 3 )*i«(TCK)»Er,EflGT/ 2 .*Ty 
i- ! 'i i'-r.! iJ +rfi2.*i)Hun*CL*Ul+POl*2.*iil^4.*Rl**n+2,*POl*C2,* 

1 ■ Mi Hi ) J 4 -t. 1 '( 2 ,»ANURM(l)=^:Rl *»2 + 2 .*ANiJRM(l))- 

2 ( I ,-U*PiiU<'AnUHM( 1) JJ*iiA(T(K)*ELENGT/2,^Ux 

521 l=-sy.:i l+E*( 2 .’f'UHUj 4 *CE=»'R 2 '- 4 .*Rl»» 2 »R 23 + 4 .*POI»ANOHH(l)*Rl»R 2 
1 (-.;,^i.’‘A(i')HfU 23 «-Plt<' 2 *Cl,- 4 .*POI)=*‘ANORM( 2 ) 3 »WTCK)*EGENGT/ 2 .*Uy 

51 U'=.'.i :f 2 + Ett 2 .«UHl>NXt(POI»K 2 « 4 ,*Rl»* 2 »R 2 )+ 2 .»POI*CANURM(l )» 

1 ( 1 +H.i t-AfJOR?l(. 2 J^Rl»’!' 2 )+E*C 2 .^ANORM(l)=^Rl*R 2 +A«ORMC 2 }))» 

2 '.TU'. J-*‘!::Le^'*GT/2.*Ux 

fl 2 .> 7 = 5212 +F’f t 2 ,i‘UKy(J*(P 0 i*Rl- 4 ,»Rl=*'R 2 ** 2 )+ 2 ,»P 01 »(AN 0 RM(l)* 

1 !i 2 -^*- 2 +/ifiURfU 2 )»iU*H 2 )+e'*‘( 2 .!>ARORM( 2 )*Rl*R 2 l-ANDftMCU ))=»= 

2 J’!M K j •I'ELbAGT/.^ .’^uy ' ' 

51,2 2 = 51 22 +F»C 2 .*;J:bJflt(E*Ri- 4 ,^Bl*R 2 =^» 2 )+ 4 .»POl»AMORMt 23 »Rl*B 2 + 

1 2 , )*R 2 » 1 ' 2 -( 1 ,"' 4 ,»P 0 I)»AN 0 RMC 1 ))’^«T(K)»ELENGT/ 2 .*UX 

52 2 /=S 2 22 l-t''*t 2 ,*i>Ki 1 *J«‘(E*H 2 + POI* 2 .*R 2 - 4 ,i 1 'R 2 *’»^i)+ 2 .*POl*C 2 .* 

1 ,uiUl!'!t 2 )*f< 2 <‘+ 23 +Ef C 2 .»AiJOFM( 2 )=i'R 24 ^‘ 2 + 2 .»ANORMC 2 ))- 

2 Cl U*PUI3*AfiOR:i(2) J»w'ftK)*EGERGT^2.4Uy 

20 Cfi'lTLnilF, 

Rr::'l !.0<IJ ; fl'iD 

if.;( rt i: f.^i:^ 3 f;f::lt 44 it 4 **$!t.^* 4 ‘*t***’t:*^^t^**^***^¥$***f^******* 

* Contrinutions trom ELSIRS are asssembled in this routine 

SUHliiiu i, i ih: STH5Lr'rb;)Ub,lMUB/lNDE,XN,ifN,XI,n,MI,TRACT,DISPL,NC) 

* XfliHilClt Rl':Al.*h<./\“-.i, 0 -Z 3 

DiMEfiSlON Xrity?) rY;n 97 ) ,XI( 20 ) |YIC 20 ) 

Dl^lENSiUM GEPt J 2 )r-^rt 32 ),IMDli( 10 ),INDEa 0 ).RC (23 
D 1 MEM 5 IUM JiSPL(l^^i)fTHACT( 193 ),BVCl 93 ),KnDE( 193 ) 

COMf.'UtM/a/GFP, ri',P| 

COMMU(^/B/GE,Pi) 1 

IB 0 D. 1 C = 1 




X y 


DO 20 Isl -NX , , „ 

iSHAU=G. ? 2 MA 12 = 0 . ?ZRA 22 = 0 . 

DU 15 J=1,MCC1) 

DO 10 K=lfNDUB 

IF(J.EQ,X?n;HCFJJGG TO 15 

D-3EG=d ; .NEMD=J + 1 
MPHE=HCC XBODY-l) 

IFXtrilJDY.Ey.l) NPREsO 

r, F t Abn i ■ ,GT . NC ( XBODY 3 JNEDD=NERl)-f]C C 1BU0Y 3 +NPRE 


M'Mi U.G,REND,XItI),'atI) fXrUflBEOfYNCNBEOy ^ , 

. : t ! ; ' . o )■ y H C n E D I) J , D 11 i , n 2 U , D 1 1 2 , D 2 1 2 , D { 2 2 » D 2 2 2 , 
.<.,1 i ,52li,5J 12, S212,Si22, 0222, TRACT, DISPL) 

V t = ^.:,ui1. + Dl 1X40211-5111-5211 


li ; * \ A- ^ 


„ ; 4 : j 1 1 2 4 D 2 1 2 - S 1 1 2 - S 2 1 2 
^ » 2 2 4 • ) 1 i 2 + D 2 2 2 - S 1 2 2 - S 2 2 2 


15 

25 


cr;:,’.:, i. 

I ■ i ' ’ ‘ i - I 


! ,2’:0 .1,a1 U) ,Y U1J ,ZmA 1 1 ,ZMA12 ,ZMA22 
( '.A, 1, ,5lJ 4.51 
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20 Cii 


f Ktr-f ,it: ttitt Doiindary are calculated below 

ur< U*i Cl j 

Q,:.,u ! ! dUYt'MMU; ,l,.'UB,Hn)ii,Np,XN,rN,TRACT,DlSPL,flC3 

* ) id.,H(jUNU rBOUM2,BOUN22 

-lou Hi! .t(. i('y,j a.hx 

49 C!!'.". il.'Ji: 

; dHD 


♦ HiUj! y routine Is used tor caiculatind stresses at the boundary. 

♦ Cal inn ro.ii.lfiH is dlCBUY 

.ad'. dJtJB,lNbB,lNDB:rND#X»Y,T,U,NC) 

» IfinljiClT K(-,mI:4*H(A-H,IJ-’Z) 

;U97 j ,Yt‘i7) ,1(193) ,IJ( 193 J ,NCC2J ,GFPC32),WT(32) 

Ul.!i:..oii)i, ItiUlKll I, IliDhClO) ■ ■ 

Cndnnl /.i/UFP,ivl,PI 

♦ Wl.'T'l* (0,790) 

V9C For: ;Vi t mx, 'HHOGKAd IN .BOUNDl SUBRUUTItiE'l 

bO J' 1 = 1, ,!il)lti:l ■ 

30 iFt Ji.Pl.d) .KO.fd I GO TU 20 



*k<j 

t>0 


iH' 

70 

10 


IF OOdiiFO OF Ooru-'ACES change following statement 

SHOULU BE CHAlOOiO ACCOHDINGLY 

NCON1=0?MCU02 = OC(1 i 
DO 20 I=NC0N1 + 1 ,f*CUN2 

IFd.NF.IinjGO TC 20 
DO 40 J=1,ND0H 
IFCINi)6Cd),EQ,0LKiU TO 50 
GO 'r'O bo 
Fsi 

i -j. 1 = 1-1 

] n;,<,.si = I + 2 

1 1 1 1 f. 3 

ij'' 1.'-.' 'h- 


iv; 

1 ■- ;, f 

U . , .'.. C J o. H ) lrtINl=NC0N2 

Jl (’ »! . . I'.'f’ / ■“ i '3 iHuUS2=NC(Ji'i l4 1 
!•' I . ■ C ! >•; H-.,'- TO 10 
I’v, •(;. ■'i + ; 

x:"'. C'u 3 ♦/ 

:>■* ; i, 

ii -ij. ,1-}. ; 

I „ > ..Pi, ns 1 - ; n 2=2»lHLuSJ 



1 ? 


c. 


t.,;,,. . ■>. .XtK) -X(XPLUSl J,X(lPLUS 2 }, 3 f{IMlNl),YCKJ, 

•I ' ^ ‘ , r ^ ^ 5 , ti ( i ^ 1 ) , li 1 1 u , 1 ] c n 1 ) , T 1 1 1 3 , u c 1 0 2 Tr K C 1 2 5 , 
^ * ■’. rn,VC.i.2J /"iD.'jij ,BDUM2#fGUW223 - 

* y J ,Bnui‘^12,BGUN22 


2U 


Ci. 

KT' 

L'-;) 


if 

: ‘I; i i . ' 


i ■f h r i-***-*^ 4 * **’^**** 4 -* ^****»^'¥***^*t* ******** **^***** 

♦ t in IS tne streKf;e5 eacfi node and passes tnein to BNDRY 


20 


tin,*.. 'nil I Li.f SlPbnY t XMlIVl , XO , X 1 . X?/ YMIN i , YO 
' f ^22 f T:: , J-.'iAl I f ZtiA 12 f ZMA22 ) 


1 ')2' ,').;i 


fYlfY 2,U10fUll,U12 fTlf 


r IT RKAli^tU A-b^n-Z) 

Dlni Cne 17, / J rWKSt?) ,CVt7) 

C; >'■ ' lb. I / 'i /Gii f Pi> 1 


t.PiVG 

hi/ 

ti. 

DC.) V.' 
CVdJ 

n ^ ”1 '"I i •> 

4 * 


:p JKT t (. At -X : Hid ) »*2 + ( XO-YMINI 3 ) 

:puKrUAi“X(') *>l‘2 + ni-Y 0)*=*2) 
:o.,iHTa X2-Xl3**2 + tY2-Yl)»»2) 


■I/} 


UM d,' ■ ! "i. t I 

Cf'F(^,J)=C. 

CfJBTXPlJt- 

C V ( I .) = U UU ~U1 0 3 / I bLENO 3 + ( U 1 2-Ul I ) / f ELEN 1 ) 3 *0 . 5 
Cvr2 J = UD21-D2<.)/CELtNU3 + CU22-U2l3/CELENU 3»0.5 
CV 0 3=11 ; CVt4)=12 

XpsC (. Yl-Vei/ELErtl + t Y0-1MIN1)/ELENO)»0,5 
XN=:t CXu-XI 3/ELtM + tXMIM-X03/ELEN0 3 *0.5 
COFd f i 3 = -Y.\;CC)t U ;23=XN 
COF ( 2 , i 3 =- YA ; C IJF O , 4 3 = XN 

C0FCi,i>3=XN;Crit (.Hbisy'i 
C0FO,6 3=X,J;(:iit (4,7 )=:YN 

sn=i./(2 0ci.+P'iij^‘Gt:J;si2=-poifsii 
COFOf U = -l.,’CfJt t5,-5 3 = ail;COt C5 
CDFia, 'i3 = -l,;COFt6,53=BU;COF(6 
CDFC/,23= 


, 7 3 =s 1 2 
,73=S13 

I = - 1 , r C U F 1 7 , 3 3 = - 1 . ; C D F ( 7 , 6 ) = I . / G E 


CA!iD FOiHXt 
C/.M, (••• 


•'0 i HXF ( 7 , CCt f 7 , wKB . UP f 0 3 
■‘'lAYFC?,! ,CGt fl ,mS,CV ,1 ,Q) 


ZA'd 1=CV(53 
ZAA12=CV t(j3 
2-'A22=C7t7 3 


i;i I' 

k ■ 

♦ ^ ^ » if ♦ '4 i! 4 




» 

* 

* 


* 


c.ii'i s-.d.ect preferred no, of p>oints of Gaussian quadrature, 
t'rn.-di 1/ or tny node ano the Internal point decides the no, ot 

I Sr'-: ij «'■' 1 i r, i ' t’- S 

F' X7 ,V,>,,m,\i,,x 2,Y2,NGP,GFP,S!»1') 

DV ; f'r, -.d"- ( .cv 3 f 



• «« n ’ * ' « ’ *•»* 
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10 


1% / I ^ . I i.V'' I I; A Z* 

H( ( C> 1 1 a 2)/2--APJ =»*2 + ( (yi + 5f2^/2,"y^)**2) 

K!. )P{:(u/.-.xut‘l2+cn-yi)»^2) 

> !,'l -i ; H L*, ij / I'iJ 

U ( i . Ii> ,>/[ . , 2 ) NUP=:2'i 
!!• t? '' L’l s.’I.'tI . , 'I J ;vbP=lb 
If i:,.M . :A,->zr2 

xrft- ..'-ij -■'.■•’sio 

IFII'-'C’ .U} 

l! .'VjPsft 

'-iGPsS 

;.UP=4 

Seipct: 1 .)ri Pas Peen restrictea to 12 Gaussian points. 

If vPriat'ie I'.aussidn points are needed , suitacle modifications 

are 'lee in., 3. 

NGP =12 

GPP(I )^'!,1252 J3 J'-aib ; ( I )=0.249l47045y 

GFP(^3si~s,FPC J ) - r«'I(2)s:WTC13 

GFP( i J = J<>7«31 4V-3PV82 ; ( 3)=0. 23349253654 

GFPt t.l = -GFP( 3 ) ;i?-a(4)=:ls(TC3) 

GP P (5) = ' ,587 3 J 7954 1 ,* aT { b ) =0 .20 3 1 67 426723 

GFFCnjs-uFPtbJ - ?aI(63=WTC5) 

GFPC / 3 = v,76y9u2b742 ? iaI ( 7 ) =0 . 1 6007 8 328543 

Gf’P(8) = -i,;FPC7 J ?fel(8)=WT(7) 

QPP {s )xi\<:ii,.\Hl2bb S ? ,*W1 (9) =0,1069 39325995 

GFF(2f Js-GFPfyJ rA'i(203=WT(93 

GFPC113=0, 9815606342467 ; Wl ( I 1 )=0 . 047 1 753363865 
GFP(12)=-GH’PU1) ?W'i’(12) = »4'nil) 

FoUO'.^lng !4AG routine gives the Gaussian points i the weiqnts. 

CALL i;>61iii3FCD0li3AZ,-l,,l,,l,NGP,WT^G6'P,aJ 


RETUK 


LUf) 


Fortnuia given uy smith ^ Llij are used tor calculating 
stresses at tne Internal points, 

SIJBJ^UUTlwe 6XAG‘l{.CUNLen,PEAKfXI,yi,Nl3 

IMPLICIT Rf:Ar<*^8CA-!i,0-2.J 

OlMEhSinh Xlt2'M,Vl(2<'!),GFP(32),WT(32) 

Cn.’''r.|U(J/A/GFP, W,l ,Pi 

.'-'515124,10) 

H.'.)r:'ATtlOX/ 'EXACT SOLN.') 

00 1 = 1,141 

yiM) = YllU-0, 086295664 

A'- 5 = )*F2 + YI(I)tYltl) 

M f'lLio^xxtx) )4'*2+yi(i)*yici) 

P fUK2/AKl) 

01 =i uC *30531 (2 , tFAC+ ( AK2+AK 1-4 .»COriL£rj ♦♦2 3 / AKl ) 

P '*=y:«.(.5,-FAC)/tAKl»DEN3 

P tPAC)/UKi»DEN) 

S' • ii \r;*U(X)/PI)»t tC(jl4LEf'i^*2 + 2,*XHlJ*XI(.I) + 

1 2 . «y r ( i ) + y.xii) j4.pH2/CGWJji£N ' ■ 

2 ■.-,Fj»i,*xif I)4PHl 3 

Ly, . t I i. i ' / t ) ^ (.CG.g4tO*P'i2-XX C I 3 *PtTl 3 
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,/ I . 


» V, I C 1 1 » M( I J 


«!l I ;i. U' i , i '.01 fiJll'.l i ,&TH1?,,STH22 

15 C‘ 'A , w Ji'.l i'! . 'i ) 

20 Cr^-''; i.ivJ'' 

f-'.i 'i -ii!) 

n< ■ • - ♦ 

» All tne roit laes tran. follow ^re for constant elements. * 


xn 


Oj calculator* stresses at the internal points. 


laiiaa C 0 1 »< a. C 0 I , X I , X I # XE , XE , TRACT ,DISPL , MC , M, HE) 

* Xr'/.>,. u'i i’ ,n';al.r»H(.;»-'l,.J..Z) 

I'l O' iOiiiO .tU2'-) , YU20) ,XEC97) ,YEC97) , TRACK l9i),DlSPL( 193) , 

COO'.'*"/ s/.:P:;,Pt!i 

CO'M.iU 

lYpTiorj , 7 90,) 

790 FOP. '.MCinX, 'PHiKOOia IN INTRNL SUBROUTINE*) 

Ire- i.!-o-o)Gn Tfi 105 
DO 1''5 t;=:i,ra 

'/iM/U. 1st' , ; ;Z./1A22=0, 

»»»l IF N" 1*H.H UK SURFACES chance: FOLLOWING STATEMENT 
•*?} SHODini si: CHANGcL 'ICCOKDINGLX 

DO r-i- 1=1. acu) 

C WRITERS, ♦)i,ii 

IF C' “U na . i*!u io:? 

302 XF( l-UCCU ) iO 3,3<U, Un 

304 TfUiPl =XiUU + l);TF'MP2 = VEtJfU 

XECuni^xuU); Y&t.i+i ) = iE(U 
GO VO -tt.’h • ■ 

303 DO l"u i:-=2f;l 

" ^ IFIU-UCUE) ) i' Of i' 7, 306 

107 TEMPlsXKCUC (K 1 f I J ; rEMP2=XtlNCCK3+l) 

‘ ■ XECUFU = xi;ut:tK-u + i ) 

XE(J + U = YEl i<HK-l J+1) 

GO TO 305 • 

306 CONTINUE 

GO ''^U 30 5 

IDl XnC..’E + lJ=XFCl); YtslUE + l ) = YE11) 

3‘>5 Call stkesscxj u) . 'n,(UfXEij),YECJ),xEij+i ) keu+i), 

1 ')! il ,0211 ,03,12 # 0212,0122 ,D222,Slli , S2 U . 3l I 2 1 S212 , 51 22 . S222 , 

2 -'III,::!:) 

DD KK=j,n 

irc.ucy.NCtKF) 3GU TU 309 

GD '!'>» 308 • 
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btiAiUib rMi,. .'iTRb.SMXi4,yN,X£l.yF.l,XK2,iF:2, 
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2 P () i , .; J 
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KOI.' /Pl(HPtH,JJ 
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SI ! 1 =(■ . ? 32 1 1 =0 , ; Sll 2=0 . ? 32 12=0 . ? SI 22=0 . ,* S222=0 . 
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00 2‘' K = i,ucp 

XWCK) = tXI.?tXbl )/2, + GFPil<)’^iXK2-XEl 1/2,0 
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Mf' '?*■,. j.iO t C2) »R1 442) fE^C^.^AWCHMC n 4HdP24A!JUKtH2) ) ) 4 

>r(. ■ ) n'.i, : 'TdR. 



21 


C 

e 

20 

♦ 

$ 


i 

!• /: 4 1 
V. ■i' ( > 

ol. 

r»2> 

AOfil' 

(i.- 

f >' i 

Cor. 
HI.l 
♦ ^ * 


i i'i 

MIm’ 

U 

UI' 


■•*>• H , . * I >ri Oi . ♦ c Vfj 1 4 R 1 , 4 HI 4S2 4 4 2 ) ■+ 2 . 4 pu 1* c ANORM ( i ) ♦ 

j- H ■■.J«‘H!'‘'!V2)+RM2,4AUGRH(2J4HJ4H2 + A!'H:iHm(1)))4 

’•'vr / ''* ' " ' ■ 

M' ♦ V *: , ( fcip J -4 , 4Rl »H2442 ) f 'i.* Pni*Ai’^r)KMC2)4R14H2t 
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GO TO 113 • 
cnriTiMUE 

GO TO 113 

XfiCRt + l ) = XE(1 ) 1 ybtiiE + U = ytU ) 
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f FC ULO.riCCKR) ICO TO 117 
GO- -M 116 ■ 

KF(J41 J=TEMPl,*yE(J4l )=TEMP2 

iVLioiE ■ ■ 
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H3 2 = Hr2«CRA'niPi‘ (201 1 *lO )-Cri# 
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Datas ■ ' '■■vl**’apr over a steel toundation. 
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*■ tniinilo lounaotion. «hen radius ot the wheel 
:i 1 L it.'its ioti !j .ire also cnanoed by the same factor, 

'-I re U'lo clo'iiental lencjts are very small when 
(..ir , a no. IPS. Jnts to define the contact region 

■' ill) ohI: ti'f , 


24 


9 '5 

9,2*1-'. 
9.4; . 
6 
1 
3 

Q 

11 

49 

51 

56 

ys 

•-y.36v' >, 

-y.29<'^' 
"•'•,24<6 . 
, 20(:< 
•’9,160' ' 
•"V , 1 2 9 m ' ; 

-0,80*6 i 

-U, 40<>', . 
6,00i'"> 
0,4 357 1 


. t 


' -• + ■ - 


If I 


K * ' 


. 1 

1 .^ 

t M f ’ f- 4 rl 

i 

■ I Mt..i 

V M ■!• <, 1 

I .'*h+/' 0 
: t M) fi'M ’ 
» ‘ i ■ 4* M I !' 

7n I 


9, 4 35''7 'iVOtL-i 
-0,4357 ;"7ri+iM 
-0.4 3577>-^7:fK' i 
-0, 16577''3 !'■»■' 2 
-«,3l2a. l'3r’ + 02 
-0.4244 4“ If'+f’V: 
-O.488ootai; + 02 
-0,498097 3n+o:; 

— 0,498t'y7J!(+fi2 
-0,498i. 97 3D+02 
-0,498i.!y73D+«!2 
— 0 . 4 4 8 7 1. 5 1 1) + 0 2 
-0,3501 8?HO+it2 
-0,2l2o Hin+u2 
-0,17101010+02 
-0,12910950+02 
-0,86824>'9D+91 
-0,435778/0+01 
-0,34878?'io + **l 
-0.19974670 K i 
-0,1799ollo + )-l. 

™9,15y77;>7M«.ri 

-0,l39-o'l /■'+* ’ 

-<1,999 M 1 7 4-i t< I' 
-0 . 799 ■o.*'. ■■'+'" 

— *.1,599"' ; **. -> I ' + < 1 ' . 
-9,399'''f6 

-»*. 199''i*' -■ i i tc.'. 
0,199 
9 . 399'-'* ■■O' 7 


0 5 9 9 ' i ' 

ol'/yy-'o.' 

U. 999 6 4 


0 

I * ' 4- ^ * 
1^. 


0* i 19 I 

0.1390.1 1 /o+;;i 
0, 159*577 /!-+"\ 
o.ny'-'.:,! H'+r-i 
0. 199<''467:' + ( i, 
0.3487’{2-iO + ot 


o V20.;. M + i-.i 9,76 567 JOO+.-iS 

'-9'>'>6(:>i.tD + 0u '} .30'.tOuOoo+01 


^33 10+91 


• : , Ti ‘-C -i 4 . . 

' .'I '’17821-)+02 
' ,15n2 3;-iO + 92 
•- ,7: :*• ^nt5!';+02 
3 i3 + 02 

*■ , 3 ',*1 6 4.1 i> + y 2 

■7. j5i'4d9ao + o2 

'■•.tOi !225o+?2 

' .5 J* 99-'00 + 92 
o,49Hi,, 97 4n + 02 
■•'.OqH,; y7 ?n + 02 
o.y')Hf’'ji740 + 0? 

.99 8.' 9740+02 
9,97171749+02 
'* .B9*'6t>920+f)2 
6.76429220+02 
u,6;jB&i58D+02 
54357790+02 
9,54357790+02 
*». 45612210+02 
0.45b4221l)+02 
0,27940150+02 
0.14310780+02 
0,47189540+01 
0,30153690+01 
0,17«^37u90+01 
0.75961230+00 
9,19026550+90 
0,121798 ii) + '*''' 
*i,39905f 90-ni 
0. 3239oO9n-0l 
<i.2559»ni'>-f(1 
H.ta598-75;)-ni 
■•'.I 4 '*"91 l')-01 
<: ,999''797o-02 
'•* 599-02 
“ '59 30-02 
9120-02 
■9620-93 
f H 2 ** 3 


1 5 ^ 

i-. f 

4 :|- 
IfJ 

1 6 ^ 


ryrii,^u2 


5 ^ 

10-7 I 

'.',1959 8750-01 
-;.2559H3 3D-01 
6 , 32 .3*'*6a9n-C 1 
*' .3?0'»5M9n-C'l 
9.12179840+90 


15 I 



25 


u,435'-' /< 

'-,171- •' } 

' , 212 . ,i . ? 

< . 35 '’’ ; .■• Hf ’- ti.v 

*' , 4*^^* '»? V' + f '; 
fi , 49 Mi '( 73 ' I ft !'/ 
V,49Hi (7’5i!+! 

C , 4 9 B 1 ,« 7 ^ • f- 1 . V 

, 4 B B I ' > ( 1 ■ ' f t ’ 2 

C ,‘* 24.1 !■' 4 tif , ;> 

^ , 1 *J i’ ■ ‘ '■ < > 0 ft ! ’} 

0. 174'^ J?yii + <.'2 

*! M.sij 

'0 , H 7 2 1'l !■' 1 1> ’ ' f * i. 

0.4507 3?f!,ff: 

5. 

0,200. .;M f.f. 1 
O.lROiMu . rjfdt 
0 , 1 6 C I til < ‘ ' f t ^ 
U, 14t'' (■''< 'i+t. I 
U, 120.;., I i'Off’l 

0, lo< ' ..oont!.: 

0,80''i.(i'!< 

0.60t ! 'M*' 

0, 4t)0iM.)( l.i)f('»! 
0, 20<i< f((J< 

- 0 , 200 .. h ' O ' 0 D +( m . 
-0, 40Onnt!t.*[)f III,. 
-0,60('Ofjt'O0+(i0 
- 0 , 80 OihUiOI>+(,)O 
»0, lOOooOOD+Ol 
-0.120*'0(0.'0+01 
-0, 14»'(ifH,»on + 01 
-0,16C>»»o0On+i)l 
-o.isoouC'OD+oi 
-0,20000000 + 11 1 
-0.36900590+01 
-0,4567 J740+O1 
-0 , 872*^o4i>f.i+( 1 


‘ 1 2 jn + jn 

’.'.17! 77i}a!)t.-a 
5‘ Ib.i^Oi'. + oi 
' . 47 i s't5 if. + Ol 
'!,:i i Jj '‘ 7 HI .+02 
' -27 9,i<. i5r,f02 

0. 456422 !r 
!'.45G'’221 n+02 

1 . ', b ‘ 1 J 577 9 nf 02 
0 .54 357790+02 
' .o.’OSlf; 8ntl.'2 
' • / 6 '4.2922r’'+'t2 

< 0920+02 
'.97171741. +02 

or.j 

<-•. f. 00 + 00 
f-f.oonf oo* 

9 0 + 00 

' mD+OO 

0,1 <>l I . ..f AO + lJC 
( ,o:.i <■ {) ( i ’ n f 0 (j 

'( ,< •' o n+OO 
*! . I 0 ( H ' l ‘<; r ) +'‘'<') 
o.'.’0'’Coi' on +00 
".ooOi.'OfAip+ou 
<‘1 (•(*( or^+oo 

i : ,..( 9 r (:/>(, tjn + f'O 
<: , ^rif.Atioon+oo 
<*,"!)"'*6{><m+0'j 

f.o^ !'OOf.on+fH) 
0 ,CO(H' OOOD + CO 
0 ,Of'00<<n(jfi+no 
(j,oci<'»of.ob+oo 
U.OO'-OOOOO + OO 
i> ,»,'0<iO(if»on+00 
O.OOOOOiiOO+OO 
O.OooooGun+OO 
O.OOOOOOOD+OO 
o,OfM»ou(jOi)+oy 
o.oocoof.ob+oo 
0 , 00000 (»on +00 
o .yo^^'ccoo+oo 
o.oo«‘0''i*nn+oc’ 
u.onco 'oon+Co 

(> ,). fii.OiHitjO + OO 


-0.174532 

li. + ! 

^ S ' 1 ’1 » ‘ 

!*iir»+o»> 


■’ 

* ; ^ i u 1 f '■ n !: i 

1 • 'J n + 0 0 

-d . 1 (ti I « . ' 

'i-t-r i 

i 

!.)r)+00 

-0. 1 .1' • ' > 


mtJ ^ y t 1 U } 

''9d + C2 

«tj. Ini' II' 

4-, ^ 


/' '»5 + 0i 

O.dO'-' 

r -ff « 

«il; ^ t ’ ’ 

l.nn + CJ 

0.1/0 , 

j 

• 1 * * 

! . p + 0 3 

'-’.li-' ! = 


*’• ^ S * " ' ’ ' ' 

’"•*,''’^ + (*2 


&oooooon+oo 
, l? 0000 !>n +02 
, I 2 unooun+o 2 

. UOOOOOOD+OO 

,ooooooun+ou 

.ooooooon+oo 

.OOOOOOOD+OO 

, 07000000+00 


9 ■ !. ('1 !)+t;i' 1 0,09000000+00 

i 9 ■/(‘. iD-ot i o . ooooooon+bo 


1 

t 

1 

t 

t 

1 

1 

t 

1 

1 



27 


7 5 

*: ^ i 


t 

o.ooooooon+f'o 

1 

7 A 


^ i * [ .) 

i 

".ooooooon+co 

1 

/ 7 

'■'5 1" 1 

^ ^ 

M * >.« 

i 

0 ,0f)0oooon+oo 

1 

/ A 


'St ,f , .j ♦ V 

1 

0 , ooonooon+oo 

1 

7q 

^ i 

1,. -1 - 

1 

o.ooooooon+no 

1 

6 

i"." * ^ 

i ^ f • . ^ 

i 

0,0000090n+90 

1 

8i 


taU. 1 i_t A ■ ) +'0.‘ * 5 

i 

o.oooooooo+no 

1. 

H 2 


^ A \ A » j +• A 

t 

o,ooooooon+oo 

1 

H 


iilA 

1 

O.COOOOQOO+OO 

1 

H 


t‘»*» i 

1 

9.ooooyoun+oo 

1 

y S' 


{ »-,..■{•/' ,}^|' /■ 

1 

9,00000000+90 

1 

8 6 


t i 0- ;t 

i 

0, 3 0000000 + ilO 

t 

87 

‘ • * 

’ 5 ; A A. -) 

1 

0,900000UD+00 

1 

8 8 

■ : . ^ 

‘■At ? A4-4 ’a' 

t 

o,ooootioon+oo 

1 



J + Aa 

1 

O.OOOOOOOD+OO 

1 

90 


. I't A < i A A 5 i ) f A^ A‘ 

1 

O.OOOOOOOn+oO 

1 

91 

1"- ^ i 


1 

U. 60000000+00 

1 

9 2 

* # < 

IM} AfA; ;; j) f *; 

S,i 

O.COOOOOOO+OO 

0 

93 

r .f 

A Af i A A';f’'- 

0 

u ,o<''0oooon+oo 

0 

9 4 

'v» • f. 

Ajf , aAi 6^ 1) fi;. 

>J 

o,ooooooor+6o 

0 

9S 


. « A ' M / j i , A IJ i t ^ ’ 

1 

»A 00000000+00 

1 


74 

30 

75 

20 

7r. 


7 7 

07 

7h 

26 , 

79 

25 

no 

2 1 


93 

;'l 

2 7 

83 

21 

P4 

20 

P5 

1<> 

at, 

IH 

'07 

17 

3o 

16 

73 

31 

72 

32 

71 

r] 

7' 

3 ' 

6 '» 

35 

6 8 

36 

6 7 

37 

h(» 


65 

30 

61 

4'. 

6 3 

41 , 

1 i ^ t 

62 

i 04 It f 

A'A 

t,i 

11 

o' 

-A 

4 L” (I 

59 

3 

+ 5 

t • . 


O.H4591nFi4( i 
0,206264HiM' t 
o.iooooof r n 
0.2291831 u-t**' 
O,11459lS0 + <:? <: 

0,206264«0+‘-l C 
O.t 0000000 + ^-2 0 

f} b 6 


I 37b.,99u+(!i 

1 0‘, < i' <.'O:}+02 

4S'U662u + 0r; 
1375.: 990 + 01 
.229i8ill) + 0l 
, l5C0t0OOU + 02 
f) 6 


0.0«754y3D+00 
0. 1604282U+01 
0.4000000U+01 
0.2000000D+02 
0.6&75493D+0U 
0.1604282D+01 
0.4000000D+01 
0. 70000000+02 


0.91673240+00 

0.18334650+01 

0.50000000+01 

0.91673240+00 

0.18134U50+01 

0.50000000+01 



Results from CO^Agt'^r* 


to 12: 


E HO, DISPL^CEHEHt X DlS|>LACEHEHf f‘ Tft&CtlOH X tJRACTION Y 


.i. - 


. . jmn . i ?>'><> 1 F4.n0 

2 

-0.58299E+00 

-0,l229lE+00 

3 

-0.57544E+00 

-0.17403E+00 

4 

-0.57544E+00 

-0,17403E+00 

5 

-0.55982E+00 

-0.21273E+00 

6 

•0.51783E+00 

-0.28638E+00 

7 

-0,44909E+00 

-0.34750E+00 

8 

-0.36069E+00 

-0,38365E+00 

9 

-0.32311E+00 

-0.38957E+00 

10 

-0.32311E+00 

-0.38957E+00 

11 

-0,27199E+00 

-0.38970E+00 

12 

-0,27!99E+00 

-0.38970E+00 

13 

-0.16589E+00 

-0.36009E+00 

14 

-0.ni671E-01 

-0.29918E+00 

15 

-0.!9ft23E-01 

-0.20922E+00 

16 

-0,79944E-02 

-0.19101E+00 

17 

0,15083E-02 

-0.15094E+00 

18 

O.88537E-02 

-0.11751E+00 

19 

0,13247E-0! 

-0.77104E-01 

20 

0.14055E-01 

-0.66935E-01 

21 

O,13fi47E-0I 

-0.44724E-01 

22 

0.13788E-01 

-0.4070lE-0i 

23 

0.13520E-01 

-0.36064E-01 

24 

0.12830E-01 

-0.30346E-01 

25 

0.10989E-0! 

-0.23008E-01 

26 

0,8991 3 E-02 

-0.16807E-01 

27 

O.73899E-07 

-0.11740E-01 

28 

0.67345E-02 

-O. 77991 E-O 2 

29 

0.66352E-02 

-0.49718E-02 

30 

0.65300E-02 

•0.3259311-02 

31 

0.63377E-O2 

-0.31844E-02 

32 

0,54625E-02 

-0.48204E-P2 

33 

0.37396E-02 

-0.75830|-02 

34 

Q.16425E-0? 

-0.11478E-0I 

35 

-0.28106E-03 

-0.16514E-01 

36 

-0.10197E-02 

-0.i9800E-0l 

37 

-0.13936E-02 

-0.22036E-01 

38 

-0.16369E-02. 

-0.23702E-01 

39 

-0.17963E-02 

-0,24985E-01 

40 

-0.18670E-02 

-0.25971E-01 

41 

-0.26797E-02 

-0.28520|-01 

42 

-0.27181E-02 

-0.27936E-0i 

43 

-0.50124E-02 

-0.165841-01 

44 

-0.88181E-02 

0.28727E-03 

45 

-0.14756E-01 

0.19179E-01 

46 

-0.22955E-01 

0.392O1E-O1 

47 

-0.73727E-0} 

0,11 i46E+00 

48 

-0.14948E+00 


49 

•0.25123E+0O 

0.19452E+00 

50 

-0.25123E+00 

0-12452E+O0 

51 

-0.30236E+00 

0.19462E+00 

52 

-0,30236E+00 

0.19462E+00 

53 

-0.34098E+00 

0.188761+00 

54 

-0.43512E+00 


55 

-0.5i38OE+00 

0.80499E-01 

56 

-0,56867E+00 

-0.18292E-01 

57 

-0.104i2E-03 

-0,90102E-04 

58 

-0.1S735E-03 

-0,34838|-03 

59 

-0.20480E-03 

-O. 7 O 869 E-O 3 

60 

-0.22345E-03 

-0.92692|-03 

61 

-0.24041E-03 

“0,11I13E-02 

62 

-0,25405E-03 

-0.14040E-02 

J 

-0.26484E-03 

-0,27971E-03 

-0,14992E-02 

-0.17754E-02 


fi 

O,O0O00E+OO 
O.OOOOOE+00 
0,O00OOE+00 
O.OOOOOE+00 
0,OOOOOE+O0 
0,00000E+00 
O,O0OOOE+0O 
0,O0000E+0O 
-O,5OOO0E-O1 
-0,50000E-0! 
0,O0000E+00 
O,OOOOOE+O0 
o,oooooE+6o 
O,OO00OE+OO 
0,00000E+00 
0,OOO0OE+OO 
O.OOOOOE+00 
0,OOOOOE+OO 
O.OOOOOE+00 
O.OOOOOE+00 
O.OOOOOE+00 
0,OOOOOE+O0 
-0.71970E+00 
-0,1 1906E+07 
-0.15794E+0? 
-0,19527E+07 
-0,l3642E+02 
-0,49990E+0J. 
0,€0004E+QQ 
O.14075E+02 
0,1 9257E+02 
0,15603E-t’0? 
0.11713E+02 
-0,34051E+00 
0,OOOOOE+O0 
O,0000OE+00 
O,00000E+00 
0,OOOOOE+O0 
O.OOOOOE+00 
0,O0O00E+00 
O.OOOOOE+00 
O,00O00E+00 
O.OOOOOE+00 
O,O0OOOE+00 
0,O0000E+0O 
O.OOOOOE+00 
O.OOOOOE+00 
0,OOOOOE4-O0 
0,OOOOOE+O0 
0,OOOOOE+O0 
0,OOOOOE+O0 
0,0O00OE+00 
0,OOOOOE+O0 
O.OOGOOE+00 
O.OOOOOE+00 
O.OOOOOE-i-OO 
0,O0O00E+00 
O.OOOOOE+00 
O.OOOOOE+00 
0,OOOOOE4-O0 
O.OOOOOE+OO 
O.OOOOOE+00 
O.OOOOOE+OO 


. .A A.AA A AP.xA A 

-0,12000E+02 

-0.12000C+02 

0,OOOOOE+O0 

O,OO0OOE+0O 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

0,OOOOOE+O0 

O.OOOOOE+00 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

0,O0O0OE+O0 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

0.53175E+00 

0.^9765E+02 

0.39485E+02 

0.48817E+02 

0.55643E+02 

0.58802E+02 

0.59940E+02- 

0,55698E+02 

0;4«l43E+02 

0.39008E+02 

0,292P3E+02 

-0.8ES29E+0a 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O,OOOO0E+O0 

0,OOOOOE+O0 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 

O.OOOOOE+OO 










Jttfpg 

, 1211 ’ 

..Hi 

' 1*1 

l*llli!i"8i 

0«2S583E'*03 

®*llfl?E-03 

9,2p04E-03 

0»3|653E»03 

0,29054E<oa3 

Q.18256E-03 

0.13f60E-03 

, ^A..JUSJ*AAAm.AA. 

•0,143S9E-03 

0»OOOOOE4>00 

©•OOOOOE+00 

0*00900E4>00 

O.I3392E*03 


342 E -03 
;-03 






if « 0 .d^. 


f 11*1^1" 
fi•2l 

*'0*193?S1«*02 

d*22422E-92 
0,2l983E«>02 

•o!2iji||2p| 

-0j24O70E-O2 

-0.22716E-02 

-O,2t240E-O2 

-9t20223E-O2 

• 0 , 194 e 0 E - O 2 

•0;i8837E*02 

- 0 ; 15546 E -02 

-0,144|7E*02 

-0*11356E-02 

•0.94377E-Oi 

•0,ail97E-03 

-0.35601E-03 

jwyfit. 

-0,12573E-03 

o,oooooe^oq 

o » oooooe 4 >oo 







:8:|S|??i:8l 

■8:88888118? 

e *|| i 43 E 42 

9.|f527E-»-02 

0.157|4E402 

0.1|904E+02 

0.2i270E+00 

0,O0000E+00 

O.$OOOOE+0O 

O.OOOOOE+OO 

o.ooeooE+oo 

O.OOOOOE+OO 

O.OOOOOE+00 

O.OOOOOE+OO 

O,O0090E-f00 

O.OOOOOE-t-00 


O.OOOOOE-fOO 

0,17978E+00 

0,«5725E-03 

-0,17580E+00 

O.OOOOOE+OO 



If 


- O . lstfSE+OO 
0 . 90 l ® eE^OO 
O . OOOOOE +00 
0 , OOOOOR + O 0 
0,00800E+00 
O.OOOOOEtOO 
O.OOOOOE+00 
O.OOOOOE-^00 
O,O®0O0E+OO 
O.OOOOOE+OO 


O.OOOOOE+00 

0,23|15E+00 

0,64858E<t>00 

0.2i984E4'00 

0,00900E+00 


Results froB tite second Pr^ram ( Page 1 3 to 221 
Stresses at internal points ) 


DDE 


X-COORD 


Y-COORD 


SIGMA-X 


SIGMA-Xy 


SIGMA-Y 


1 

2 

3 


6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 
19 
2 C 


O.PCfOE+Ol 
■0.509|®i+C 
0 , 50000 E + « _ 
•Cl500C0E+0i 
0^50000E+01 
•OlsOClOOE+Oi 
O.SOOOOE+Ol 
-O.5O0|OE+Oi 
0.5Wwi+Ol 
-O.SOWE+Ol 
O.SOWE+Oi 
O.OCOOOE+OO 
O.oOOOOE+00 
0.OO9fOE+OO 
O.OOOOOE+00 
OloOO ' OOE + Of ) 
0 ld 0000 f !+00 
O . OOOOOE^OO 
0,oOOO<)E+Of» 


O,2Oe00E+01 
*I§?DOE+Ol 


_ OOOE+Ol 
0.40000E+C1 
0,60000E+01 
O.PtOOE+01 
d,80dOOE+Ol 
0.90000E+Oi 
0 * 100091^02 
O*l§0Ooit62 
0.1|00©E402 

8'liSSSI+?i 

OlioOOOE +01 

0*60000E+Ol 

O*8OOO0E+O1 

O,t000OE+02 

0,15©OOE+02 

0.20000E+02 

0,25OOOE+C2 


.65255E+01 

,55S23E+01 

0.51275E+01 

0,51853E+0l 

0,32030EtOl 

0,34l89E+0l 

0,17830E+01 

0,20408E+01 

0*94216E+00 

0,11S79E+Ol 

0la|l47E-0l 

O*l||66E+00 

0 : 4 i 225 E+Cl 

0*8||74E+Oo 

0 : ltS 65 E - 0 l 

*|190§E400 

!S1102E+00 

a,45636E+00 

0.51336E+00 

0.54148E+00 


8 


0*32421E-I-C1 
0.23836E+01 
0.49685E+01 
■0-42778E+01 
0,45483E+01 
0,43065E+01 
0,36509E+01 
■o’ 56575E+01 
0.28484E+01 
0*29555E«i-01 
l5858Et0l 
l7208Et01 
O,26615E+0l 
0*l2670E+01 
67285E+00 


V” IP 

4 : 


,40857E+OO 
-0*27272E+00 
-0.I2782E+0O 
-0*73942E-01 
-0,48353E-01 


:8 


:8 


0.i6065E+f 
0,i0020E-l-0. 
O.493i0E+Ol 
.37295E+01 
. esIssE+oi 
0,54614E+01 
0.684|8E+0 
,60522E+0 
,65848E+| 

8:I8HII:8 
:8:ll«|lt8 

"oltlllal+l 
"ollJfilito 
:E:ltl!li:Si 

-0*40505E+01 


